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(57) The invention relates to a method of processing 
a digital signal comprising samples (a^ a^ a$) each 
represented by a coordinate on each axis of a coordi- 
nate system with n dimensions, where n is the dimen- 
sion of said digital signal, characterized in that it com- 
prises the following steps: 

determining among the coordinates of said samples 
the extreme coordinates (L, R, B, T) on each coor- 
dinate axis (x and y) of the coordinate system, 
for all the possible combinations between said sam- 
ples taken in pairs, determining a minimum distance 
among ail the distances calculated along each axis 
between, on the one hand, the coordinate of a first 
sample, and, on the other hand, the coordinate of 
each of theelements ofthe group comprising a sec- 
ond sample and copies of this second sample ob- 
tained by translation of said second sample in dif- 
ferent directions and according to a value which de- 
pends on the extreme coordinates already deter- 
mined. 
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Description 

[0001] The present Invention relates to a method of 
processing a digital signal comprising samples each 
represented by a coordinate on each axis of a coordi- 
nate system with n dimensions, where n is the dimen- 
sion of said digital signal. 

[0002] The invention applies in particular to the field 
of images. 

[0003] The invention applies generally to the case of 
a digital signal which is not necessarily an image and 
which comprises a set of samples representing physical 
quantities. 

[0004] Thus, the digital signals concerned by the in- 
vention may for example be video signals, audio signals, 
(music, speech, etc.), computer signals, signals output 
by a facsimile machine, multidimensional signals (multi- 
band imagery), etc. 

[0005] When it is desired to encode samples of a dig- 
ital signal, and more particularly of an image, the con- 
ventional techniques of encoding are essentially based 
on the approach which consists in taking the samples 
forming the digital image in a predefined order, for ex- 
ample line by line and from left to right, and of encoding 
the amplitude of these samples. 
[0006] The samples of the image signal concerned by 
the encoding may be pixels of said image signal. 
[0007] Nevertheless, methods of decomposition of 
digital image signals are in widespread use today of the 
type where the decomposition is into frequency sub- 
bands. 

[0008] The value of these methods derives from the 
aim of concentrating the maximum of information Into a 
predetermined sub-band, in this case the sub-band of 
lowest frequency. This sub-band then contains coeffi- 
cients which constitute a subset of the samples of the 
digital image, this subset preserving the essential char- 
acteristics of the digital image. 
[0009] Thus the digital images may be, for example, 
decomposed into sub-bands and the amplitude of the 
coefficients of each of the sub-bands is then encoded, 
[0010] In this case, the processing order of the digital 
samples is considered as implicit. 
[0011] It should be noted that it may be chosen to en- 
code only certain digital samples among all the pixels of 
the image or among all the coefficients of the sub-band 
of the decomposed image. 

[0012] Furthermore, it may be desired to adopt an op- 
posite approach to that already described on encoding 
samples of a digital signal, for example of an image. 
[0013] More particularly: considering that the ampli- 
tude of the digital samples is implicit, I.e. that it conforms 
more or less to a predetermined model, a description of 
a path between the digital samples concerned of the sig- 
nal will be sought, this path then being encoded. 
[001 4] In order for the encoding to be carried out with 
a minimum throughput, it is necessary for the descrip- 
tion of the path between these samples to be as com- 



pact as possible. 

[0015] If the digital samples of a bi-dimensional signal 
are represented for example by points in a Cartesian co- 
ordinate system comprising two axes of coordinates, it 
5 is then necessary for the path linking ail these points to 
be the shortest possible to minimize the number of bits 
necessary for encoding. 

[0016] In practice, the description of this minimum 
path involves the calculation of a distance between all 

10 the points taken in pairs, the definition of this distance 
being based on a vector linking two points. 
[0017] In order for the path to be described between 
these points to be minimum, known solutions of the 
"traveling salesman problem" are relied on, this being a 

15 well-known problem in mathematics. According to this 
problem, a traveling salesman must pass through a cer- 
tain number of towns in an order which minimizes the 
total distance. 

[0018] In the case of interest here, solutions to this 
20 problem provide a minimum path joining the points or 
digital samples considered. 

[0019] However, to solve this problem, the concept 
used of distance between two points must first be de- 
fined. 

25 [0020] If M 1 and M 2 are two points with respective co- 
ordinates (x 1( y,) and (x 2 . y 2 ) in a Cartesian coordinate 
system, and V is the vector of displacement between M 1 
and M 2 having the following coordinates: 

30 

vx=x 2 - Xj 



vy=y 2 -yi> 

35 

the distance d between the points M-j and M 2 , also 
called norm, Is usually defined by one of the following 
three types: 



w I — 5 2 

d = Vvx + vy 

or 

45 

d=max(lvxl,lvyl) 

or, 

50 

d=lvxl + lvyl. 

[0021] When the distance has been selected, the so- 
lutions mentioned above to the "traveling salesman 
55 problem" provide a list of vectors V between the different 
points to be passed through along the chosen minimum 
path. 

[0022] With the chosen vectors having a norm of low 
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value, It may be understood that the encoding of the cho- 
sen path is advantageous in terms of throughput. 
[0023] Nevertheless, the applicant has realized that it 
would be worthwhile to find a new method and a new 
device making it possible to further reduce the number 5 
of bits necessary for the encoding of samples of a digital 
image. 

[0024] On solving this problem, the applicant has 
found, in a surprising manner, a new method and a new 
device for processing a digital signal which provide, for 
each pair of samples, the smallest possible distance be- 
tween those samples. 

[0025] This new processing method and device serve 
more particularly as a basis for the creation of a new 
method and device making it possible to reduce the 
number of bits to be encoded in the signal. 
[0026] However, in a more general manner, the new 
processing method and device found by the applicant 
may be used independently for other applications. 
[0027] The object of the present invention Is thus a 
method of processing a digital signal comprising sam- 
ples a 1t a^ a 9 each represented by a coordinate on 
each axis of a coordinate system with n dimensions, 
where n is the dimension of said digital signal, charac- 
terized in that it comprises the following steps: 

determining among the coordinates of said samples 
the extreme coordinates L, R, B, T on each coordi- 
nate axis x, y of the coordinate system, 
- for all the possible combinations between said sam- 
ples taken in pairs, determining a minimum distance 
among all the distances calculated along each axis 
between, on the one hand, the coordinate of a first 
sample, and, on the other hand, the coordinate of 
each of the elements of the group comprising a sec- 
ond sample and copies of this second sample ob- 
tained by translation of said second sample in dif- 
ferent directions and according to a value which de- 
pends on the extreme coordinates already deter- 
mined. 

[0028] In a complementary manner, the invention re- 
lates to a device for processing a digital signal compris- 
ing samples a 1( a 2 , a 9 each represented by a coor- 
dinate on each axis of a coordinate system with n di- 
mensions, where n is the dimension of said digital sig- 
nal, characterized in that It comprises: 

means for determining among the coordinates of 
said samples the extreme coordinates L, R, B, T on 
each coordinate axis x, y of the coordinate system, 
for all the possible combinations between said sam- 
ples taken in pairs, means for determining a mini- 
mum distance among all the distances calculated 
along each axis between, on the one hand, the co- 
ordinate of a first sample, and, on the other hand, 
the coordinate of each of the elements of the group 
comprising a second sample and copies of this sec- 
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ond sample obtained by translation of said second 
sample in different directions and according to a val- 
ue which depends on the extreme coordinates al- 
ready determined. 

[0029] Thus, in other words a "frame" containing the 
samples considered is defined in the coordinate system, 
"copies" of this frame on each side of it are made by 
translation, adjacent to said frame, in different directions 
and the smallest possible distance is determined be- 
tween a first sample of the frame and a second sample 
of said frame, as well as between this first sample and 
the "copies" of the second sample inside the respective 
copies of the frame. 

[0030] In this way a new displacement vector W is de- 
fined for each pair of samples so processed and of which 
the norm is the smallest there is. The norm of the vectors 
W is In fact less than or equal to that of the displacement 
vectors V of the type described above. 
[0031 ] According to one feature, the method compris- 
es a step of determining a minimum path between the 
samples a v a 4 , ag, a 7 , a 5 , e^, a 3 , 89 from the different 
minimum distances determined previously along each 
axis and for all the possible combinations of pairs of 
samples. 

[0032] Having so determined new displacement vec- 
tors for each pair of samples considered, it is ensured 
that the minimum path obtained by applying a known 
solution of the "traveling salesman problem" to these 
new vectors will be less than or equal to that which would 
be obtained with the above-mentioned vectors V using 
the same solution. 

[0033] A minimum path thus shortened results in vec- 
tors of reduced norm which are thus more advanta- 
geous to encode in terms of throughput. 
[0034] It will be noted that the new displacement vec- 
tors W may be determined before the application of the 
above-mentioned solution of the "traveling salesman 
problem" or simultaneously with the implementation of 
that solution. 

[0035] According to one feature, the value which de- 
pends on the extreme coordinates corresponds sub- 
stantially to the difference L-R, R-L, B-T, T-B between 
the maximum coordinate L, B and the minimum coordi- 
nate R, T of the samples on each axis. 
[0036] The translation is thus adapted to the dimen- 
sions of the set containing the samples in order to make 
the copies of samples outside a virtual "frame" delimiting 
that set. 

[0037] According to one feature, the copies of the sec- 
ond sample are made in the two possible directions of 
each axis and of each diagonal. 
[0038] According to one feature, the displacement 
vectors obtained from the consecutive samples of the 
minimum path are encoded. 

[0039] These vectors are next transmitted with the or- 
der in which they appear in the minimum path in order 
to be able to position the different samples after decod- 
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ing. 

[0040] The invention also relates to a method of 
processing an encoded digital signal comprising a step 
of receiving this encoded digital signal, characterized in 
that said encoded digital signal comprises in encoded 
form: 

a coordinate of a so-called initial one of the samples 
of said non-encoded signal which was obtained 
along each axis of a coordinate system of n dimen- 
sions, where n is the dimension of said signal, 
the extreme coordinates L, R, B, T on each axis of 
the samples of said non-encoded signal, 
displacement vectors W which were each obtained 
before encoding of said signal from consecutive 
samples of a minimum path determined between 
the samples of the signal, as well as the order of 
these vectors in the path. 

[0041] In a complementary manner, the invention re- 
lates to a device for processing an encoded digital signal 
comprising means for receiving this encoded digital sig- 
nal, characterized in that said encoded digital signal 
comprises in encoded form: 

a coordinate of a so-called initial one of the samples 
of said non-encoded signal which was obtained 
along each axis of a coordinate system of n dimen- 
sions, where n is the dimension of said signal, 
. the extreme coordinates L, R, B, T on each axis of 
the samples of said non-encoded signal, 
displacement vectors W which were each obtained 
before encoding of said signal from consecutive 
samples of a minimum path determined between 
the samples of the signal, as well as the order of 
these vectors in the path. 

[0042] The digital signal received was encoded ad- 
vantageously with respect to the throughput which 
makes the reception of the data so encoded more rapid. 
[0043] After decoding of the data, the coordinates of 
each of the samples of the minimum path, along each 
axis, are reconstituted as a function of the new displace- 
ment vectors, of the extreme coordinates and of the co- 
ordinate of the initial sample. 

[0044] This thus enables the image signal to be re- 
constructed in a manner just as reliable as if all the am- 
plitudes of the samples had been transmitted or else as 
if more information had been transmitted by transmitting 
the displacement vectors V mentioned above. 
[0045] According to one feature, the method compris- 
es a step of determining the position of the coordinates 
reconstituted along each axis of each of the samples of 
the path with respect to the extreme coordinates of the 
samples. 

[0046] The test carried out on the position of the re- 
constituted coordinates makes it possible to know 
whether that position is situated within the virtual "frame" 



enclosing the original samples or outside it. 
[0047] According to one feature, when this position is 
outside the frame, i.e. when the determination step re- 
veals that the position of a coordinate of a sample along 
5 an axis is situated beyond the extreme coordinates on 
that axis, said method comprises a step of adjusting the 
position of said coordinate. 

[0048] More particularly, the adjusting step consists 
of carrying out a translation of the position of the coor- 
10 dinate along the axis considered and according to a val- 
ue which depends on the extreme coordinates on said 
axis and on the position of the coordinate with respect 
to said extreme coordinates. 

[0049] As previously stated, methods of decomposi- 
te tion are in widespread use today of the type where the 
decomposition is into frequency sub-bands. 
[0050] When digital images are thus decomposed, 
the amplitude of the coefficients of each of the sub- 
bands is encoded. 
20 [0051] Given the above, the Applicant's attention 
turned to a new method and a new device for processing 
a digital signal in which the amplitude of the samples 
would not be encoded. 

[0052] The present invention thus also relates to a 
25 method of encoding a digital signal composed of a set 
of samples, with each sample there being associated 
an amplitude and at least one coordinate in a coordinate 
system, the method being characterized in that it com- 
prises the steps of: 

30 

classifying at least some of the samples according 
to a predetermined criterion, 
obtaining a curve representing the variation of the 
classified samples according to the predetermined 
35 criterion, 

encoding the data representing the curve obtained, 
encoding information representing at least one co- 
ordinate of at least one coefficient in an order ob- 
tained from the classification. 

40 

[0053] In a complementary manner, the invention also 
relates to a method of transmitting a digital signal which 
has been encoded according to the earlier encoding 
method, characterized in that the encoded data repre- 
ss senting the curve obtained and the encoded information 
representing at least one coordinate are transmitted by 
means of a network. 

[0054] The invention also concerns a method of de- 
coding a digital signal comprising coefficients, the meth- 
so od being characterized in that it comprises the steps of: 

- decoding information representing a curve of vari- 
ation of the amplitude of the classified coefficients, 
decoding information representing at least one co- 
55 ordinate of the coefficients in a predetermined order 
and associating with the coefficients an amplitude 
obtained from the information representing the de- 
coded curve and the order of the coefficients. 
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[0055] Thus, the amplitudes of the coefficients do not 
need to be encoded, so reducing the quantity of infor- 
mation to be transmitted. 

[0056] According to a preferred embodiment the dig- 
ital signal is divided into sub-parts beforehand and the 
classification of at least some of the samples is carried 
out among the samples of a sub-part. 
[0057] This division of the digital signal Into sub-parts 
makes It possible to process only some of the samples, 
and more particularly neighboring samples, thus making 
It possible to obtain curves representing the variation of 
the amplitudes of the samples having only little variation. 
[0058] This division of the digital signal into sub-parts 
also enables only small q uantities of samples to be proc- 
essed at a time, and so reduces the size of the volatile 
memory of the processing device as well as its complex- 
ity. 

[0059] In a complementary manner, the invention re- 
lates to a device for encoding a digital signal composed 
of a set of samples, with each sample there being as- 
sociated an amplitude and at least one coordinate in a 
coordinate system, the device being characterized in 
that it comprises: 

means of classifying at least some of the samples 
according to a predetermined criterion, 
means of obtaining a curve representing the varia- 
tion in the samples classified according to the pre- 
determined criterion, 

- means of encoding data representing the curve ob- 
tained, 

means of encoding information representing at 
least one coordinate of at least one coefficient in an 
order obtained from the classification. 

[0060] The invention also relates to device for trans- 
mitting a digital signal which has been encoded accord- 
ing to the encoding method, characterized in that the 
encoded data representing the curve obtained and the 
encoded information representing at least one coordi- 
nate are transmitted via a network. 
[0061] The invention also concerns a device for de- 
coding a digital signal comprising coefficients, the de- 
vice being characterized in that it comprises: 

- means of decoding information representing a 
curve of variation of the amplitude of the coefficients 
classified, 

- means of decoding information representing at 
least one coordinate of the coefficients in a prede- 
termined order and of associating with the coeffi- 
cients an amplitude obtained from information rep- 
resenting the decoded curve and the order of the 
coefficients. 

[0062] in relation to another aspect, the JPEG2000 
standard, whose description is available via the Internet 
network at the address http^/www Jpeg.org , currently 
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constitutes the benchmark method for encoding and de- 
coding a digital image. 

[0063] This method give satisfactory results In partic- 
ular from the point of view of compression performanc- 

5 es, I.e. the ratio between the quality of the restored im- 
age after decoding and the size of the compressed file. 
[0064] However, these results can be Improved on. 
[0065] The object of the present invention is to provide 
a method and a device for digital signal encoding which 

10 enable compact encoding of the digital signal. 

[0066] To that end, the invention provides a method 
of encoding a set of data representing physical quanti- 
ties, the set of data comprising coefficients, 

characterized in that it comprises the steps of: 

15 

quantizing the coefficients using a predetermined 
set of quantization symbols, 
counting the number of coefficients associated with 
each of the quantization symbols, 
20 - encoding the number of coefficients associated with 
each of the quantization symbols, 
determining the location of each of the coefficients 
associated with each of the quantization symbols, 
and 

25 - encoding the locations of the coefficients. 

[0067] The invention makes it possible to obtain good 
compression performances, i.e. a high ratio between the 
quality of the restored signal after decoding and the size 
30 of the compressed file. 

[0068] Thus, for the same quality of image, the com- 
pressed file obtained according to the invention is small- 
er than that obtained by applying the JPEG2000 meth- 
od. Reciprocally, for the same size of file, the image re- 
35 stored after decoding according to the invention has a 
better quality than that obtained by applying the 
JPEG2000 method. 

[0069] According to a preferred feature, the locations 
of the coefficients are determined by: 

40 

- the coordinates of an initial coefficient in the set of 
data, 

displacement vectors for the other coefficients, the 
displacement vector of a coefficient being calculat- 
es ed from another coefficient associated with the 
same symbol. 

[0070] The locations are thereby compactly encoded. 
[0071] According to a preferred feature, the encoding 
so of the locations of the coefficients comprises an entropic 
encoding of the displacement vectors. 
[0072] This type of encoding is simple to implement 
and gives good performance, 

[0073] Alternatively, according to a preferred feature, 
55 the locations of the coefficients are represented by their 
coordinates In the set of data. 

[0074] The locations are then determined more sim- 
ply than according to the first possibility. 
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[0075] According to a preferred feature, the encoding 
of the locations of the coefficients comprises an entropic 
encoding of the coordinates of the coefficients. 
[0076] Here too, this type of encoding proves to have 
good performance. 

[0077] According to a preferred feature, the set of data 
is a signal of the frequency sub-band resulting from a 
decomposition into sub-bands of an initial signal. 
[0078] The decomposition of a signal into frequency 
sub-bands does not In itself create any compression , but 
enables the signal to be decorrelated in order to elimi- 
nate the redundancy prior to the actual compression. 
The sub-bands are thus encoded more efficiently than 
the original signal. 

[0079] The invention also relates to a method of de- 
coding a set of data representing physical quantities en- 
coded by the encoding method already presented, 
characterized in that it comprises the steps of: 

decoding the number of co efficients associated with 
each of the quantization symbols, 
decoding the location of each of the coefficients as- 
sociated with each of the quantization symbols, 
setting each of the coefficients to the value of the 
quantization symbol which corresponds to It. 

[0080] The decoding method enables the image to be 
restored. 

[0081 ] According to a preferred feature, the decoding 
of the location of a coefficient comprises the decoding 
of the displacement vector which corresponds to it. 
[0082] Alternatively, according to a preferred feature, 
the decoding of the location of a coefficient comprises 
the decoding of its coordinates. 
[0083] The decoding of the location of a coefficient of 
course depends on the manner in which this location 
was encoded. 

[0084] In a complementary manner, the invention re- 
lates to a device for encoding a set of data representing 
physical quantities, the set of data comprising coeffi- 
cients, 

characterized in that it comprises: 

means for quantizing the coefficients using a pre- 
determined set of quantization symbols, 
means for counting the number of coefficients as- 
sociated with each of the quantization symbols, 
means for encoding the number of coefficients as- 
sociated with each of the quantization symbols, 
means for determining the location of each of the 
coefficients associated with each of the quantiza- 
tion symbols, and 

means for encoding the locations of the coefficients. 

[0085] The encoding device comprises means for im- 
plementing the above features. 
[0086] The invention also relates to a device for de- 
coding a set of data representing physical quantities en- 
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coded by the encoding method already presented, 
characterized in that it comprises: 

means for decoding the number of coefficients as- 
5 sociated with each of the quantization symbols, 

means for decoding the location of each of the co- 
efficients associated with each of the quantization 
symbols, 

means for setting each of the coefficients to the val- 
10 ue of the quantization symbol which corresponds to 
it. 

[0087] The decoding device comprises means for im- 
plementing the above features. 
15 [0088] The encoding device andthe decoding method 
and device have similar advantages to those already set 
out. 

[0089] The invention also relates to a digital appara- 
tus including the devices according to the invention, or 

20 means for Implementing the methods according to the 
invention. This digital apparatus is for example a digital 
camera, a digital camcorder, a scanner, a printer, a pho- 
tocopier, or a fax machine. The advantages of the digital 
apparatus are identical to those already set out. 

25 [0090] According to another aspect, the invention also 
relates to: 

an information storage means which can be read 
by a computer or a microprocessor containing code 
30 instructions for a computer program for executing 
the steps of the methods according to the invention, 
and 

an partially or entirely removable Information stor- 
age means which can be read by a computer or a 
35 microprocessor comprising code instructions for a 
computer program for executing the steps of the 
methods according to the invention. 

[0091 ] According to still another aspect, the invention 
40 relates to a computer program which can be directly 
loadable into a programmable device, comprising in- 
structions or portions of code for implementing the steps 
of the methods according to the invention, when said 
computer program is executed on a programmable de- 
45 vice. 

[0092] Since the features and advantages relating to 
the information storage means and to the computer pro- 
gram are the same as those set out above relating to 
the methods according to the invention, they will not be 
50 repeated here. 

[0093] Other features and advantages of the present 
invention will emerge more clearly from a reading of the 
following description, given with reference to the accom- 
panying drawings, in which: 

55 

Figure 1 represents samples of a bi-dimensional 
digital signal and their respective locations in a Car- 
tesian coordinate system; 
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Figure 2 is an example embodiment of a system 
showing, on the one hand, the processing, encod- 
ing and transmission of a digital signal of which the 
samples are shown In Figure 1, and, on the other 
hand, the decoding and processing of the digital sig- 
nal decoded according to the invention, enabling 
the samples of the signal of Figure 1 to be reconsti- 
tuted; 

Figure 3 is a view of a virtual frame C delimiting the 
samples of Figure 1; 

Figure 4 is a diagram of a displacement vector of 
known V type between two samples; 
Figure 5 is a diagrammatic view representing cop- 
ies of the frame C of Figure 3 on each side of it, in 
the eight possible directions; 
Figure 6 shows a new displacement vector obtained 
between a sample M 1 and a copy M' 2 of a sample 
M 2 among the samples of Figure 1 ; 
Figure 7 and 8 represent the minimum path ob- 
tained and the associated displacement vectors be- 
tween the samples of Figure 1 respectively with and 
without the invention; 

Figure 9 is an embodiment of a programmable de- 
vice implementing the invention; 
Figure 1 0 is an algorithm for processing a digital sig- 
nal to be encoded according to the invention, imple- 
mented in the device of Figure 9; 
Figure 1 1 is an algorithm for processing a digital sig- 
nal encoded according to the invention, implement- 
ed in the device of Figure 9; 
Figure 12 is an example of an embodiment of an 
encoding-decoding system according to the inven- 
tion; 

Figure 1 3 is an encoding algorithm according to the 
invention implemented in the device of Figure 12; 
Figure 14 is a decoding algorithm according to the 
invention implemented in the device of Figure 12; 
Figure 1 5 is a circuit for decomposition into frequen- 
cy sub-bands included in the device of Figure 12; 
Figure 1 6a is a digital image to be encoded by the 
encoding device according to the present 
Invention ; 

Figure 1 6b is an image decomposed into sub-bands 

by the circuit of Figure 15; 

Figure 1 7 is the image decomposed into sub-bands 

by the circuit of Figure 15, and partitioned into 

blocks; 

Figure 1 8 is an algorithm for approximation calcu- 
lation of the variance of the coefficients classified; 
Figure 19 shows an encoding device according to 
the invention and a corresponding decoding device; 
Figure 20 is an embodiment of an encoding method 
according to the invention; 
Figure 21 is an embodiment of an decoding method 
according to the invention. 

[0094] As shown in Figure 1 , a bi-dimensional digital 
signal, such as an image signal, comprises several sam- 



ples a 1 , a 2 , a 3) a 4 , a 5 . a 6 , a 7 , a 8 , a 9 which are for example 
here pixels represented in a Cartesian coordinate sys- 
tem with x and y coordinate axes. 
[0095] Generally, an image signal is a series of digital 
5 words, for example bytes. Each byte value represents 
a pixel of the image, here with 256 levels of gray or else 
in color. 

[0096] The size of the image or of a subset of it rep- 
resented in Figure 1 is 21 x21 pixels. 

10 [0097] The pixels may thus have coordinates which 
range from the coordinates (0,0) at bottom left to the 
coordinates (20, 20) at top right. 
[0098] The different samples form a set denoted G 
and are represented in cross-hatched squares, respec- 

f5 tively situated at the following locations: a 1 = (1 7, 2), a 2 
= (1, 3), 83 = (3, 3), 84 = (16, 6), a 5 = (1, 7), fife = (16, 9), 
a 7 = (17, 11), ag = (18, 13), et ag = (2, 20). 
[0099] It should be noted that the samples may alter- 
natively be sub-band coefficients coming from one or 

20 more decompositions into frequency sub-bands of the 
image signal. 

[0100] Such decomposition methods, for example of 
DWT (Discrete Wavelet Transform) type are widely 
known and used these days. 

25 [0101] Instead of encoding the set of all the samples 
with their amplitudes and their coordinates, the present 
invention provides for defining, via new displacement 
vectors, a series of the positions of the samples of the 
set G and to encode it for the purposes of transmission, 

30 as will be seen below. 

[0102] It will be noted that the digital signals con- 
cerned by the invention are not necessarily bi-dimen- 
sional since the invention applies for example to mono- 
dimensional signals such as audio signals or to tri-di- 

35 mensional signals such as video signals. 

[01 03] With reference to Figure 2, an embodiment of 
a processing device according to the invention is adapt- 
ed process and to encode a digital signal with the object 
of compressing it. The processing device is integrated 

40 into an apparatus, which is for example a digital camera, 
a digital camcorder, a scanner, a printer, a photocopier, 
a fax machine, a system of database management, or 
a computer. 

[01 04] The device according to the invention compris- 
es es a signal source 30, which is an image signal IM here 
comprising the set G of samples of Figure 1 . Generally, 
the source of the signal either contains the digital signal, 
and comprises for example a memory, a hard disk or a 
CD-ROM, or converts an analog signal into digital sig- 
so nal, and is for example an analog video camera associ- 
ated with an anaiog to digital converter. The image 
source 30 generates a series of digital samples repre- 
senting an image IM. 

[0105] An output of the signal source 30 is connected 
55 to a circuit 32 for determining the coordinates L, R, B, T 
of the extremes of the signal on each axis of the coor- 
dinates x, y of the coordinate system. 
[0106] Among the coordinates of all the samples of 
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the set G, the extreme coordinates L, R, B, T are thus 
determined by seeking the coordinates of the minimum 
and the maximum that the signal reaches on each axis, 
i.e.: 

L = min (x) R = max (x) B = min (y) T = 
max (y). (x, y) G G (x, y) e G (x f y) e G 
(x, y) G G 

[0107] For the set G, the extreme coordinates are : 
L= 1, R= 18,B = 2,T=20. 

[0108] As shown in Figure 3, the determination of the 
extreme coordinates enables a virtual frame C to be de- 
fined which delimits as closely as possible the set G of 
the samples. 

[0109] Returning to Figure 2, the output of circuit 32 
is connected to a circuit 34 which will provide new dis- 
placement vectors W according to the invention. 
[0110] To do this, the concept is defined beforehand 
of the distance to use for the calculation of the distance 
between two samples and, for example, the distance will 
be taken as distance d = max (lx 2 - Xfl, ly 2 - YtO between 
two points M 1 and M 2 of respective coordinates x^ y 1 
and x 2 . y 2 . 

[0111] The determination of these new vectors re- 
quires all the samples of the set G to be processed in 
pairs. 

[0112] In the prior art, it was sought to calculate the 
distance between the two samples of the pair (M 1t M 2 ) 
in course of processing (Figure 4). 
[0113] This distance associates a positive real 
number with the displacement vector V linking M 1 and 
M 2 and is also called the norm of vector V. 
[0114] According to the invention, "copies" of the 
frame C of Figure 3 will firstly be considered which were 
made on each side of said frame in the eight possible 
directions, as shown in Figure 5. 
[0115] In this manner, "copies" of the samples of the 
set G are made within the different "copies" of the frame. 
[01 1 6] The "copies" of samples are obtained by trans- 
lation of a sample considered in the two possible direc- 
tions for each coordinate axis, along the diagonals and 
according to a value which depends on the extreme co- 
ordinates determined by circuit 32. 
[0117] For example, this translation value corre- 
sponds, for the x-axis, to the difference between the 
maximum coordinate R and the minimum coordinate L 
on that axis, which gives for the two possible directions 
the respective values R-L and L-R. 
[0118] Similarly, on the y-axis the values T - B and B 
- T are found for the two directions. 
[0119] Next, the invention provides for calculating, for 
a pair of samples in course of processing, such as that 
of Figure 4, all the distances between M 1 and each of 
the elements of the group comprising M 2 and the "cop- 
ies" of M 2 , i.e. M' 2 , M" 2? M ,n 2 , M 2 < 4 >, M 2 < 5 >, M 2 < 6 >, M 2 <7) 
and M 2 ( 8 ) (Figure 5). 

[0120] The smallest distance between all the distanc- 
es calculated is then determined and the displacement 
vector W of coordinates wx, wy is thus obtained for the 



pair M-|, M 2 . 

[0121] it will be noted that the coordinates of the vec- 
tor W are written: 

5 wx = min (lx 2 - x^, lx 2 + L - R - x 1 l, lx 2 - L + R - x n l), 
wy = min (ly 2 -y^ ly 2 + B « T - y^, ly 2 - B + T - y,!), 

10 

where x lf y 1 and x 2 , y 2 are respectively the coor- 
dinates of M 1 and M 2 . 

[01 22] It will be noted that the coordinates of the vec- 
tor W supplied above make it possible to search for the 

15 sample to which the vector points or a copy of that sam- 
ple among one of the nine "frames" of Figure 5. 
[0123] Thus, for the samples M v M 2> the new dis- 
placement vector W is represented in Figure 6 and con- 
nects M 1 to a "copy" M' 2 of the sample M 2 . The "copy" 

20 M' 2 Is closer to M, than M 2 and all the other "copies" 
M* 2 , M- 2l M 2 W, M 2 ( 5 >, M 2 W. M 2 < 7 >, M 2 <*>. 
[01 24] The vector W corresponds to a virtual displace- 
ment towards a new sample M' 2 which belongs to a 
"copy" of the signal, obtained by an offset (translation) 

25 of R - L to the right (axis for which x Is positive). 

[0125] In this manner, circuit 34 determines, for each 
pair of samples of the set G, a new displacement vector 
W. 

[0126] These new displacement vectors could be 
30 used independently of the embodiment whose descrip- 
tion will follow. 

[0127] For example, when a mouse is used to move 
a cursor on a screen of large dimensions of a computer 
or other electronic apparatus, it is necessary to move 
35 the mouse over relatively long distances, which is not 
always easy given the support on which rests the 
mouse. 

[0128] The new vectors W thus make it possible, for 
example, to reduce the displacements of the mouse by 
40 defining a faster path between the control buttons of the 
mouse on which the user must click. 
[0129] It would thus be possible, by displacing the 
mouse in a single direction, to make the cursor leave 
the right side of the screen and to make it return virtually 
45 immediately on the left without needing to make it cross 
the whole screen and thus without needing to displace 
the mouse as a consequence. 
[0130] Returning to Figure 2, the output of circuit 34 
is connected to a circuit 36 which will apply to the mini- 
50 mum distances calculated previously one of the known 
solutions of the "traveling salesman problem". 
[0131] Solutions of this problem may be found in the 
article of R.E. Burkard et AL, "Well-solvable special cas- 
es of the TSP: a survey", SPEZIAL FORSCHUNGSBE- 
55 REICH F003, Bericht n° 52 f Dec. 1995. 

[0132] An example of a method is also given at the 
Internet address http://itp.nat.uni-magdeburg.de/ 
~mertens/TSP/node2.html. 
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[0133] This method Is called "method by insertion" 
and relies on the principle idea which is to first construct 
a path comprising a subset of the samples to be visited. 
[0134] For example, three samples are chosen at ran- 
dom and connected in the only manner possible. 
[0135] Next, for each sample not yet inserted, a cal- 
culation is made of which pairs of successive samples 
it will best be inserted between (i.e. by lengthening the 
path as little as possible ), and this represents the cost 
of that sample. 

[0136] The sample of minimum cost is then inserted 
in the path. This procedure is reiterated until all the sam- 
ples have been inserted. 
[0137] A path of short length is thus obtained. 
[0138] It should be noted that the path is not neces- 
sarily the shortest distance possible, but the solution ob- 
tained is nevertheless satisfactory. 
[01 39] It should be noted that no solution exists at the 
current time guaranteeing a calculation of minimum path 
when the number of samples is high, and for a reason- 
able calculation time. 

[01 40] Circuit 36 will thus make It possible to obtain a 
minimum path or route between the samples of the set 
G which is represented by a series of displacement vec- 
tors of new type and which are each, by definition, of 
reduced norm. 

[0141] The vectors W |f i = 1 to 8, are the following: (- 
1,4), (0,3), (1,2), (1,2), (1,-6).. (0,-4), (2,0), (-1, -2). 
[0142] A last vector (-1, 1) serves indicate the first 
sample of the path already identified and it is thus not 
necessary to take it into consideration during encoding. 
[0143] The minimum path is written 
a-j 3436373335328339 and its length is 26 for the distance 
used, whereas the minimum path obtained without the 
invention, with displacement vectors of V type, is written 
a 1 a3a 2 a 5 a 9 a 8 a 7 a 6 a 4 a 1 and Its length is 60. 
[0144] Figures 7 and 8 illustrate by the intermediary 
of displacement vectors the minimum paths obtained re- 
spectively with the invention and with the vectors of V 
type. 

[0145] The invention thus enables a virtual path or 
route to be created In some way which a description will 
not be attempted but which will enable, by its compact- 
ness, to reduce the number of bits necessary for the en- 
coding of the signal. 

[01 46] Although the determination of the short vectors 
of new W type (determination of a minimum distance 
between two samples and between a sample and "cop- 
ies" of the other sample) has been presented as being 
independent of the operations carried out by circuit 36, 
it is possible to have a single circuit combining circuits 
34 and 36. 

[0147] In this single circuit, the determination of the 
vectors W is made simultaneously or virtually simulta- 
neously with the solution of the "traveling salesman 
problem". 

[0148] Returning to Figure 2, the output of circuit 36 
is connected to an entropic encoding circuit 38. 



[0149] Beforehand, an Initial sample of the set G with 
its coordinates is selected, for example a 1t then is en- 
coded by circuit 38. 

[0150] This initial encoded sample will be transmitted 
5 to the processing device of the encoded digital signal 
and will serve as a starting point for determining the lo- 
cations of the other samples of the set G in the coordi- 
nate system. 

[0151] Circuit 38 also carries out the entropic encod- 
10 ing: 

- of the extreme coordinates (L, R, B, T) of the set of 
the samples G and 

of the eight displacement vectors W), 1 = 1 a 8, de- 
15 termined by the circuits 34 and 36 and which de- 
scribe the minimum path a^aQajaQa^a^ag in an 
order established from the initial sample. 

[01 52] It is possible to directly encode the coordinates 

20 of each vector. 

[0153] However, if it is desired to reduce the memory 
size of the data, it is preferable to encode the list of the 
vectors, for example using a Huffman encoding or an 
arithmetic encoding. 

25 [0154] The encoding may also comprises a prior 
quantization step. 

[01 55] The new vectors of reduced norms will thus be 
encoded in a more economic manner than before. 
[0156] The encoding circuit is connected to a circuit 
30 40 in order to transmit to it a compressed file containing 
the signal of the encoded image, and, more particularly, 
in encoded form: 

the x and y coordinates of the Initial sample a 1 , 
35 . the extreme coordinates (L, R, B, T) and 

the vectors W ( , where i = 1 to 8, with the order in 
which they appear in the path starting from the initial 
sample. 

40 [01 57] The circuit 40 stores and/or transmits the com- 
pressed file to a device for processing the encoded dig- 
ital signal. 

[0158] This processing device receives the encoded 
data then using them performs the reconstruction of the 
45 compressed image. It should be noted that if an entropic 
encoding and/or a decomposition into sub-bands have 
been carried out at the encoder, the decoding device will 
perform the inverse operations to those of the encoding 
circuit 38. 

so [0159] After reception of the encoded data constitut- 
ing the encoded signal, they undergo a decoding oper- 
ation in a decoding circuit 42. 

[0160] The processing device then carries out, in a 
general manner, the reconstitution of the locations or po- 
55 sitions of the different samples of the set G as a function 
of the extreme coordinates of said set, of the vectors W| 
and of the coordinates of the initial sample. 
[0161] More particularly, the processing device com- 
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prises a circuit 44 which draws up the list of the displace- 
ment vectors W| as a function of the order of these vec- 
tors in the minimum path starting from the initial sample 
a v 

[01 62] The output of circuit 44 is connected to a circuit 
46 for reconstitution of the coordinates of the different 
samples of the set G as a function of the extreme coor- 
dinates L, R, B t T, of the list of vectors W, and of the 
coordinates of the initial sample a^ 
[01 63] The details of the operations carried out by this 
circuit will be dealt with later. 

[01 64] The circuit 46 is connected to a circuit 48 using 
the decoded data, for example to view the positions of 
the samples of the set G (decoded image). 
[0165] With reference to Figure 9, an example of a 
programmable data processing apparatus 1 00 imple- 
menting the Invention is described. This apparatus is 
adapted to process a digital signal and in particular to 
encode it. 

[01 66] Furthermore, this apparatus is also adapted to 
process an encoded digital signal and in particular to 
decode it. 

[0167] According to the chosen embodiment shown 
in Figure 9, the apparatus implementing the invention is 
for example a microcomputer 1 00 connected to different 
peripherals, for example a digital camera 1 01 (or a scan- 
ner, or any means of acquiring or storing an image) con- 
nected to a graphics card and supplying data to be en- 
coded. 

[0168] The apparatus 100 comprises a communica- 
tion bus 102 to which there are connected: 

- a central processing unit 1 03 denoted CPU (micro- 
processor), 

- a read-only memory 104, comprising a program 
"Progr 1", for the processing and in particular the 
encoding of data and a program w Progr 2" for the 
processing and in particular the decoding of encod- 
ed data, 

a random access memory 1 06, comprising registers 
adapted to record variables modified during the ex- 
ecution of the aforementioned programs, 
a screen 1 08 for either displaying the data to be en- 
coded or decoded as the case may be or serving 
as an interface with the user, who can parameterize 
certain modes of encoding or decoding, using a 
keyboard 1 1 0 or any other means, such as a mouse, 

- . a hard disk 112, _ 

a floppy drive 1 1 4 adapted to receive a diskette 1 1 6, 
an interface 1 1 8 for communicating with a commu- 
nication network 1 20 able to transmit encoded dig- 
ital data or to receive encoded data which are to be 
decoded by the device, 

an input/output card 122 connected to a micro- 
phone 124 (the data to be processed according to 
the invention in that case constitute an audio sig- 
nal). 



[0169] The communication bus affords communica- 
tion between the different elements included in the mi- 
crocomputer 1 00 or connected to it. The representation 
of the bus is non-iimiting and, in particular, the central 

5 processing unit is able to communicate instructions to 
any element of the microcomputer 100 directly or by 
means of another element of the microcomputer 100. 
[01 70] The programs denoted "Progr 1 " and "Progr 2 M 
enabling the programmable apparatus to implement the 

10 invention, may be stored for example in read-only mem- 
ory 1 04 (denoted ROM in the drawing), as shown in Fig- 
ure 3. According to a first variant, the diskette 116, like 
the hard disk 1 1 2, may contain data that are encoded or 
to be encoded as well as the code of the invention which, 

15 once read by the apparatus 100, will be stored on the 
hard disk 112. In a second variant, the programs may 
be received in order to be stored in an identical fashion 
to that described previously via the communication net- 
work 120. 

20 [0171] The diskettes can be replaced by any informa- 
tion carrier such as, for example, a CD-ROM or a mem- 
ory card, in general terms, an information storage 
means which can be read by a computer or microproc- 
essor, integrated or not into the apparatus, and which 

25 may possibly be removable, stores one or more pro- 
grams implementing the processing method according 
to the invention 

[0172] More generally, the programs may be loaded 
Into one of the storage means of the apparatus 1 00 be- 
30 fore being executed. 

[0173] The central processing unit 103 will execute 
the instructions relating to the implementation of the in- 
vention, which are stored in the read only memory 104 
or in the other storage means. On powering up, the pro- 
as grams which are stored in a non-volatile memory, for ex- 
ample the ROM 104, are transferred into the random 
access memory RAM 106, which will then contain the 
executable code of the invention, as well as registers for 
storing the variables necessary for implementing the in- 
40 vention. 

[0174] The above-mentioned data processing appa- 
ratus 1 00 may also, in a non-limiting manner, be a digital 
camera, a digital moving picture camera, a camcorder, 
a television set, a video recorder, a fax machine, a serv- 
es er, a printer, etc. 

[0175] It should be noted that the data processing ap- 
paratus comprising a device for processing a digital sig- 
nal that is to be encoded or which is encoded according 
to the invention may also be a programmed apparatus. 
so [0176] This apparatus then contains the code of the 
computer program or programs for example fixed in an 
application specific integrated circuit (ASIC). 
[0177] Figure 10 shows an algorithm comprising dif- 
ferent instructions or portions of code corresponding to 
55 steps of the method of processing a digital signal ac- 
cording to the invention. 

[0178] The computer program denoted "Progr 1", 
which is based on this algorithm, is stored in the read 
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only memory 1 04 of Figure 9 and, on initialization of the 
system, is transferred into the random access memory 
106. 

[0179] It is then executed by the centra! processing 
unit 1 03, which thus makes it possible to implement the 
method according to the invention in the device of Figure 

9. 

[0180] As shown in Figure 1 0, a first step of the algo- 
rithm, denoted S1, consists of acquiring the points or 
samples a 1t a 2 , ... ag making up the set G and of which 
the location (positions represented in Figures 1 and 2) 
is to be encoded. 

[0181] The algorithm of Figure 10 comprises a next 
step denoted S2 which performs a calculation of the ex- 
treme coordinates L, R, B and T among the coordinates 
of the samples of the set G, on each axis (x-axis and y- 
axls) of the coordinate system. 
[0182] This calculation was described in detail in the 
description of circuit 32 of Figure 2. 
[0183] During this step S2, an encoding of the ex- 
treme coordinates thus determined is also carried out. 
[0184] The following step denoted S3 concerns the 
transmission of the extreme coordinates L, R, B,T which 
have been encoded. 

[0185] It should be noted that, according to a variant, 
the extreme coordinates thus determined may be trans- 
mitted, once encoded, with other encoded data. 
[0186] The algorithm of Figure 10 comprises a next 
step denoted S4 during which the choice of a so-called 
initial location is made, that is to say an initial sample of 
the set G and of its position represented by the coordi- 
nates in the x, y coordinate system. 
[0187] This initial sample which was taken above 
equal to a lP Is the sample starting from which , at the 
device containing the decoder, the other locations or po- 
sitions of the different samples of the set G will be re- 
constructed, by virtue of the invention. 
[0188] During this step S4, an encoding of this initial 
sample is also carried out. 

[0189] The following step S5 provides for the trans- 
mission of the coordinates of the encoded Initial sample 
to the device containing the decoder. 
[0190] Here too, this transmission step may be 
grouped together with the step of transmitting other da- 
ta. 

[0191] Similarly, the encoding operations planned for 
steps S2 and S4 may also be combined with the encod- 
ing of other data. 

[0192] During the following step denoted S6, the de- 
termination is carried out of new short vectors according 
to the invention by processing the samples of the set G 
in pairs. 

[0193] More precisely, for each pair of samples con- 
sidered, the distance between these two samples is cal- 
culated as well as the distance between one of these 
samples, called the first sample, and each of the "cop- 
ies" of the other sample called the second sample. 
[0194] The copies of this second sample are obtained 



by translation of the second sample in different direc- 
tions, for examnie in the two possible directions of each 
axis and of each diagonal and according to a value 
which depends on the extreme coordinates determined 
5 at step S2. 

[01 95] More precisely, fo r each axis this last value cor- 
responds substantially to the difference between the 
maximum coordinate and the minimum coordinate of 
the samples on the axis considered. 
[0196] Thus, on the x-axis, two translations are pro- 
ceeded with of a value equal to L - R and R - L for the 
two directions concerned, and, along the y-axis, two 
translations of a value corresponding to B - T and T - B 
in the two directions considered. 
[0197] Once all the foregoing distances have been 
calculated, the least distance among them is deter- 
mined and to this a new vector of type W with a reduced 
norm can thus be attributed. 

[0198] This step is that which corresponds to the op- 
erations carried out by circuit 34 of Figure 2. 
[0199] The algorithm of Figure 10 also comprises a 
step denoted S7 during which a minimum path is deter- 
mined in known manner between the samples a v a^ ... 
a 9 of the set G, using the short vectors already deter- 
mined, that is to say the different minimum distances 
between the coordinates of the samples considered, 
along each axis and for ail the possible combinations of 
pairs of samples. 

[0200] As stated earlier in course of the operations 
carried out by circuit 36 of Figure 2, the minimum path 
between the samples is determined by using one of the 
known solutions of the traveling salesman mathematical 
problem. 

[0201] Step S7 gives rise to the vectors of type W, 
where i = 1 to 8 and whose coordinates were given 
above. 

[0202] These new vectors are shown in Figure 7. 
[0203] During the following step denoted S8, the con- 
struction is can-led out of the list of the vectors W f in the 
order in which they appear in the minimum path deter- 
mined at step S7, and starting with the initial sample 
chosen previously. 

[0204] This list of vectors W drawn up starting with the 
initial vector a t wili be used at the decoder to reconsti- 
tute the different coordinates of the samples of the set 
G and thus position these samples at their respective 
locations. 

[0205] The following step S9 corresponds to an en- 
troplc encoding of the vectors W, of the list already 
drawn up at step S8 f as performed by the encoding cir- 
cuit 38 of Figure 2. 

[0206] Follow step S1 0 provides for the transmission 
of the vectors so encoded in the order in which they ap- 
pear in the list drawn up at step S8, in order for them to 
be able to be appropriately used at the device containing 
the decoder. 

[0207] This step ends the algorithm of the method of 
processing a digital signal according to the invention. 
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[0208] Figure 11 shows an algorithm comprising dif- 
ferent instructions or portions of code corresponding to 
steps of the method of processing the encoded digital 
signal according to the invention. 
[0209] The computer program denoted "Progr 2", 
which is based on this algorithm, is stored in the read 
only memory 104 of Figure 9 and, on Initialization of the 
system, is transferred into the random access memory 
106. 

[0210] It is then executed by the central processing 
unit 1 03, which thus makes it possible to implement the 
method according to the invention in the device of Figure 
9. 

[0211] As shown in Figure 11 , the algorithm compris- 
es a first step, denoted S10, of receiving the encoded 
data making up the digital signal encoded according to 
the algorithm of Figure 10. 

[0212] During a following step denoted S21 , a decod- 
ing is made of the extreme coordinates L, R, B, T which 
were determined and encoded at step S2 of the algo- 
rithm of Figure 10. 

[0213] During the following step S22, a decoding is 
made of the coordinates of the initial sample a 1 selected 
at step S4 of the algorithm of Figure 10. 
[0214] During the following step denoted S23, an en- 
tropic decoding is carried out of the vectors Wj deter- 
mined at step S6 of the algorithm of Figure 1 0. 
[0215] During the following step denoted S24, the list 
is drawn up of the vectors W 8 in the order in which these 
vectors appear in the minimum path starting with the in- 
itial sample a 1( by virtue of the data transmitted by the 
device for processing the digital signal to be encoded. 
[0216] During the following step denoted S25, the pa- 
rameter i representing the vector W } concerned is set to 
1. 

[0217] During step S26, the calculation is made of the 
coordinates of the following sample based on, on the 
one hand, the preceding sample which is, here, the ini- 
tial sample, and, on the other hand, the first vector W 1 
of the list drawn up at step S24. 
[0218] The coordinates of the vector W, considered 
will be denoted wx. wy and the coordinates of the pre- 
ceding sample will be denoted x 1( y v 
[021 9] The coordinates x 2 , y 2 of the following sample 
are written: 

X 2 = X-j + wx, 



y 2 = y 1 + wy. 

[0220] During the following step denoted S27, a test 
is performed on the position of the reconstituted coordi- 
nate x 2 along the x-axis of the following sample with re- 
spect to the extreme coordinates L and R. 
[0221] Depending on whether x 2 is less than L or 
greater than R, this means that the situation shown In 



Figure 6 has occurred in which the following reconstitut- 
ed sample falls outside of the frame C shown in Figure 3. 
[0222] In this case, step S27 is followed by a step S28 
during which the position is adjusted of the reconstituted 

5 coordinate along the x-axis, In an appropriate manner 
according to whether x 2 is less than L or greater than R. 
[0223] Otherwise, when x 2 is both greater than L and 
less than R, this means that the sample whose position 
is reconstituted is actually within the frame C mentioned 

10 above. 

[0224] The two cases which have just been described 
are summarized by the formulae below: 

If x 1 + wx < L, x 2 = x 1 + wx + R - L 
If x 1 + wx > R, x 2 = x 1 + wx - R + L 



^ If R > x 1 + wx > L, x 2 = x 1 + wx 

[0225] It is thereby understood that the adjustment 
step consists of performing a translation of the position 

25 of the coordinate along the x-axis according to a value 
which depends on the extreme coordinates on that axis. 
[0226] More particularly, the value depending on the 
extreme coordinates corresponds to the difference R-L 
or L-R according to the position of the coordinate of the 

30 sample considered with respect to the extreme coordi- 
nates. 

[0227] During afollowing step S29,atest is performed 
on the position of the reconstituted coordinate y 2 along 
the y-axis of the following sample with respect to the 

35 extreme coordinates T and B. 

[0228] When the test carried out shows that the coor- 
dinate y 2 determined at step S26 is greater than T or 
less than B, then step S29 is followed by a step S30, 
during which an adjustment of the position of that recon- 

40 stituted coordinate along the y-axis is made, in a manner 
adapted to the position of that coordinate with respect 
to the extreme coordinates T and B. 
[0229] In fact, the results of the test mentioned above 
signifies that the reconstituted coordinate of the f ollow- 

45 ing sample is located outside the frame C illustrated in 
Figure 3. 

[0230] Otherwise, when the reconstituted coordinate 
y 2 is less than T and greater than B, then the sample is 
one of the samples of the set G and not one of the "cop- 
so ies" of these samples which were made uniquely to im- 
prove the encoding in terms of througput 
[0231] The tests made on the coordinate y 2 to deter- 
mine its position with respect to that of the extreme co- 
ordinates and the operations carried out during the ad- 
55 justment step are summarized below: 

if y n + wy > T, y 2 = y t + wy + B - T 
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lfy 1 + wy<B,y 2 =y 1 + wy-B + T 



If T > y 1 + wy > B, y 2 = y 1 + wx. 

[0232] In a similar manner to that stated above con- 
cerning the adjustment step S28, according to the posi- 
tion of the coordinate y 2 with respect to the extreme co- 
ordinates B and T, a translation of that coordinate y 2 re- 
constituted at step S26 is carried out with a value B - T 
orT - B according to the case envisaged. 
[0233] The algorithm of Figure 11 next comprises a 
step S31 during which a test is carried out on the value 
of the parameter i with respect to the value N identifying 
the last vector of the list of vectors W,. 
[0234] When the result of this test is negative, step 
S31 is followed by a step S32 during which the param- 
eter i is incremented by one unit and step S26 already 
described is returned to. 

[0235] Otherwise, when the result of the test carried 
out at step S31 is positive, a step S33 is proceeded to 
during which the results obtained during the preceding 
steps are brought together, so enabling the set G to be 
reconstituted as shown in Figure 2. 
[0236] This step ends the algorithm of Figure 11 and 
the method of processing a digital signal according to 
the invention. 

[0237] With reference to Figure 12, an embodiment 
of an encoding device according to the invention is 
adapted to encode a digital signal with the object of com- 
pressing it. The encoding device is integrated into an 
apparatus, which is for example a digital camera, a dig- 
ital camcorder, a scanner, a printer, a photocopier, a fax 
machine, a system of database management, or a com- 
puter. 

[0238] The device according to the invention compris- 
es a signal source 50, in this case image signal IM which 
is either a fixed image or a sequence of images. Gen- 
erally, the source of the signal either contains the digital 
signal, and comprises for example a memory, a hard 
disk or a CD-ROM, or converts an analog signal into dig- 
ital signal, and is for example an analog video camera 
associated with an analog to digital converter. The Im- 
age source 50 generates a series of digital samples rep- 
resenting an image IM. The image signal IM is a series 
of digital words, for example bytes. Each byte value rep- 
resents a pixel of the image IM, here with 256 levels of 
gray or in color. 

[0239] An output of the image source 50 is connected 
to a circuit 51 for classifying digital words or coefficients 
obtained by their decomposition into sub-bands by the 
circuit 51 . 

[0240] The circuit 51 enables a classification to be 
made according to a predefined increasing or decreas- 
ing order, the order not being, of course, absolute. It will 
be seen In what follows that certain compromises may 
exist in the classification. The operation of the circuit 51 



will be detailed below with the help of the flow diagram 
of Figure 13. 

[0241] The circuit 51 is connected to an encoding cir- 
cuit 52, which will determine, as a function of the sam- 
ples, a curve representing variations of the digital words 
orcoefficients. The encoding will also associate with this 
curve information that represents at least one coordi- 
nate of the coefficients or digital words in the image. The 
encoding of the classified coordinates may also com- 
prise a quantization and an entropic encoding such as 
an arithmetic encoding or a Huffman encoding. The 
curve as well as the classified coordinates form a set of 
encoded information. 

[0242] The encoding circuit 52 is connected to a cir- 
cuit 53 for processing of encoded data for transmitting 
to it the compressed file CF containing the encoded im- 
age. The circuit 53 stores and/or transmits the com- 
pressed file CF to a decoding device 54. 
[0243] The decoding device 54 receives the encoded 
data then performs the reconstruction of the com- 
pressed image from the curve and the classified coor- 
dinates. It should be noted that if an entropic encoding 
and/or a decomposition into sub-bands have been car- 
ried out at the encoder, the decoding device will perform 
the inverse operations to those of the encoding circuit 
52. The decoding circuit 54 is connected to a circuit 55 
using the decoded data, for example to view a decoded 
image IM'. 

[0244] With reference to Figure 13, the encoding 
method according to the invention will now be de- 
scribed. 

[0245] The algorithm as described in Figure 13 is for 
example stored in the memory 1 04 of the device 1 00 of 
Figure 9 in order to be subsequently loaded into the 
memory 1 06 for its execution. The storage medium stor- 
ing the algorithm is integrated or can be detachably 
mountable on th edevice 1 00. For example, the storage 
medium is a floppy disk or a CD-ROM. 
[0246] The digital image IM is decomposed at step 
301 using preferably, but not limited to this, a decompo- 
sition into sub-bands as described with reference to Fig- 
ures 15, 16a and 16b. 

[0247] The coefficients thus obtained are next, at step 
305, divided into blocks. This division into blocks having 
been performed, the central processing unit 103 of the 
device 1 00 of Figure 9 will, at step 31 0, take a first block 
and process it independently of the other blocks. The 
central processing unit will take the blocks for example 
in the following manner: the first being the block situated 
at the upper left of Figure 17, the second being that 
which is at the top to the right of the first and so on. 
[0248] In variant form, the processing by block will be 
carried out sub-band by sub-band following an order 
similar to that already described. 
[0249] This operation achieved, at step 315 the cen- 
tral processing unit will carry out the calculation of ap- 
proximation of the coefficients contained in the block by 
a function. This calculation divides up Into different sub- 
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parts described with reference to Figure 18. 
[0250] A first step 801 of Figure 1 B will consist in an 
encoding of the sign of the coefficients. The coefficients 
being of real value, the sign information will first of all be 
encoded independently, by transmitting for example a 
binary value at 0 for a plus sign and at 1 for a minus sign. 
[0251] Other methods of sign encoding are possible 
such as those mentioned in the JPEG 2000 standard 
mentioned above. 

[0252] This operation achieved, and the coefficients 
then being considered as absolute values and thus all 
being positive, step 805 will consist in classifying the co- 
efficients of the block from largest to smallest. 
[0253] In a variant form, a threshold value S1 could 
be fixed, below which the coefficients are considered as 
not significant. Putting in order will then be carried out 
as far as the last coefficient of value greater or equal to 
S1 , the coefficients considered as not significant are not 
processed. 

[0254] It should be noted that other variants are pos- 
sible, for example from the smallest to the largest. 
[0255] In what follows it will be assumed that P coef- 
ficients are classified {that is to say that P coefficients 
have a value greater than or equal to S1). 
[0256] This operation having been achieved, the cen- 
tral processing unit 103 goes on to the step 810 which 
will choose, from the series of P coefficients arranged 
from the greatest to the least and named by the series 
x tl x 2 , x 3 ,... ,x P , a set A comprising a finite number of R 
parameters a t , ^ a R . For example, a decreasing ex- 
ponential family could be chosen, written f A (x)=a v exp 
(a 3 .x+a 4 )a2. 

[0257] At steps 815, 820, 825 and 830 the method will 
determine among a family of real functions that which 
reproduces as accurately as possible the evolution of 
the values of the x, as a function of their rank /. 
[0258] At these steps, the parameters a., , a 2 , a 3 , and 
a 4 will be determined such that the function f A is the clos- 
est possible to the points (1 , xj, (2, x 2 ),...,(P, x P ). This 
typical problem is called function regression. Numerous 
techniques exist for fin ding the best set of parameters. 
For example minimization of the root mean square error 
RMSE could be sought, equal to 



[0259] To achieve this, it could be chosen to quantize 
each parameter a, (step 830) into a finite number of pos- 
sible values. There are then a finite number of different 
values of parameters of A. The RMSE value is then cal- 
culated at the step 820 for each set of parameters and 
is then stored with the value of the parameters associ- 
ated at step 825. When all the parameters have been 
incremented (step 81 5), the central processing unit will, 
at step 835, choose the set which minimizes RMSE from 



among the parameters stored at step 815 and stocked 
in the memory 106. 

[0260] Other faster methods are described in D. Mo- 
hda et AL, "Minimum complexity regression estimation 
5 with weakly dependent observations", IEEE Transac- 
tions on Information Theory, Nov. 96. 
[0261] Let Aopt be the name of the set of parameters 
which minimize RMSE. 

[0262] This operation achieved, the central process- 
10 ing unit then goes on to the step 320 of Figure 13 and 
will encode the parameters of Aopt in order to store them 
in a file representing the encoded Image, or transmit 
them to an addressee via the communication network 
120. 

15 [0263] At step 325, the central processing unit will cal- 
culate the path which minimizes the cost of encoding, 
that is to say find the path linking Q elements of the set 
{x 1( ...,x P } such that the cost of encoding of this path is 
minimum. A measurement of encoding cost is thus 

20 made available. 

[0264] In typical manner, this measurement is ex- 
pressed, for a given signal, In the form B+X.D, where B 
is the rate associated with the signal, D is the distortion 
with respect to the original after decoding of the signal, 

25 and X is a parameter fixed by the encoder or the user of 
the encoder. 

[0265] The parameter X determines the balance be- 
tween compression and quality. The distortion is a sim- 
ple measurement to calculate. This is because, when 
30 the decoder receives the information about the localiza- 
tion of the k 4 * 1 coefficient, the decoding value will be f A (k). 
[0266] Thus, at the encoder, when the k th sample is 
in course of being encoded, the distortion is measured 
between the value of this sample and f A (k). 
35 [0267] The rate is also simple to measure since we 
have an entropic code which is associated with each dis- 
placement vector. The number of bits necessary for the 
encoding of a vector can thus be obtained directly. 
[0268] This type of problem is well known in combina- 
40 torics. It is known as the traveling salesman problem, in 
which a certain number of towns have to be gone 
through in an order which minimizes the total distance. 
In this case, our problem is a particular instance called 
the time-dependent traveling salesman problem. 
45 [0269] Indeed, the cost associated with the passage 
from one sample to another depends on the path al- 
ready followed, since the decoding distortion of a sam- 
ple depends on the rank at which it will be decoded. 
[0270] The solution of this type of problem is found by 
so use of an evolutionary algorithm. An example of such a 
method can be found in X. Wang et AL, "The traveling 
salesman problem with separation requirements: an ex- 
amination of alternative formulations", or else in R. Burk- 
ard et AL, "Well-solvable special cases of the TSP : a 
55 survey", these two references being available at the In- 
ternet address http://citeseer.nj.nec.com/cs. 
[0271] Once the optimum path in terms of rate-distor- 
tion has been obtained, the central processing unit goes 
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on to step 330 which consists of encoding the first co- 
efficient. 

[0272] it should be noted that it is necessary to en- 
code the first coefficient transmitted in a different man- 
ner to the others. In order to simplify the encoding of this 
coefficient, one could for example arbitrarily encode its 
location using a vector of displacement with respect to 
the point of coordinates (1 6, 1 6), that is to say the center 
of the block. 

[0273] Next, for all the other coefficients except the 
last, the central processing unit will encode the displace- 
ment vector at step 335. 

[0274] At this step, there are transmitted the coordi- 
nates of the vector of displacement from the coordinates 
of the previously encoded coefficient. 
[0275] Displacement vector is understood to mean 
the vector of which the coordinates are the differences 
between the coordinates of the second point and of the 
first. It should be noted that in theory this solution may 
seem less beneficial than a conventional solution in 
which the position of the coefficients would be directly 
encoded since a difference over 32 coefficients is en- 
coded over 6 bits instead of 5 (from -31 to +32). Never- 
theless, the nature of the signal is such that in general, 
the coefficients of high amplitude are found to be accu- 
mulated in the same place in a sub-band. It is therefore 
worthwhile to encode the points with reference to each 
other, such that the displacement vectors remain very 
small on average. 

[0276] Thus, by applying an entropic encoding to 
these vectors, it will be possible to obtain a code that is 
more compact. 

[0277] For as long as all the vectors have not been 
encoded, steps 335 to 345 will be reiterated. 
[0278] Once all the coefficients have been encoded, 
the central processing unit goes on to the following step 
350, which consists in verifying is all the blocks have 
been processed. 

[0279] If they have not, at step 355 the central 
processing unit will take the following block from the 
RAM memory and reiterate the steps 31 0 to 355. 
[0280] If they have, the central processing unit has 
terminated the encoding method, and the encoded in- 
formation is then stored or transmitted via the commu- 
nication network 120. 

[0281] In variant form, the encoding of the coefficients 
may be carried out in the following manner. By taking 
for example a block of size 32x32 coefficients, all that is 
required is 5 bits for encoding each coordinate, in rela- 
tion to the x and y axes. 

[0282] From a parameter Q less than or equal to P, 
the coordinates of the x, are transmitted or stored in the 
order x 1 ,x 2 ,x 3 ,... until x Q . 

[0283] By modifying the value of the parameter Q, it 
would be possible to vary the compression rate and the 
quality of the image to be compressed. The more Q is 
small, the more the image will be compressed, but the 
higher will be the degradation of the reconstituted signal. 



All the blocks of all the sub-bands are encoded in this 
manner. 

[0284] According to another variant form, instead of 
encoding the coordinates over 5 bits it is also possible 
5 to use an entropic encoding such as a Huffman encod- 
ing or an arithmetic encoding. 

[0285] According to a further variant, the encoding is 
identical to that already set out in steps 301 to 355 ex- 
cept that no variation of the rate distortion parameter is 
implemented. 

[0286] According to another variant form, the sign en- 
coding is not necessary and for example the encoding 
as already described will only be carried out on the pos- 
itive coefficients then independently on the absolute val- 
ue of the negative coefficients. 

[0287] With reference to Figure 14, the method of de- 
coding conforming to the method of encoding previously 
set out will now be described. 

[0288] The decoding algorithm is stored in the non- 
volatile memory 1 04 or the hard disk 1 1 2 of the device 
1 00 of Figure 9 and is loaded by the central processing 
unit 1 03 into the random access memory 1 06 of the 
same device. The storage medium storing the algorithm 
is integrated or can be detachably mountable on th ede- 
vice 100. For example, the storage medium is a floppy 
disk or a CD-ROM. 

[0289] The algorithm comprises 10 steps numbered 
401 to 445. The decoding device constituted by the de- 
vice 1 00 containing the decoding algorithm receives the 
information to be decoded via the communication net- 
work 120, or will read the encoded information from one 
of its storage means. 

[0290] At step 401 , the central processing unit will 
read the first block to be decoded. It should be noted 
that the encoder and the decoder use the same block 
division. 

[0291] in the same manner as for the encoding, the 
blocks will be processed independently from each other 
in a sequential manner, in the same order as for the en- 
coding, but it may also be carried out in a parallel man- 
ner. 

[0292] At step 405 the central processing unit 1 03 will 
read the parameters of the approximation function as- 
sociated with the block to be processed and decode 
them if necessary. 

[0293] It should be noted that the decoder uses the 
same family of functions as the encoder 
[0294] The following step 410 will consist of reading 
and decoding the first sample of the block, then reading 
the displacement vector if the encoding has encoded the 
information in vector form. 

[0295] Once this operation has been accomplished, 
the central processing unit will attribute, to the coordi- 
nate coefficient obtained by the displacement vector, the 
value of the amplitude given by the position of the coef- 
ficient by virtue of the approximation function reconsti- 
tuted from Its parameters. The decoding of the sign In- 
formation will also be carried out so as to be able to at- 
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tribute its sign to each coefficient decoded. 
[0296] This is performed at step 420. It should be not- 
ed that if the encoding had been carried out without dis- 
placement vectors but in a simpler manner by encoding 
the position of the coefficient in the block, the calculation 
of the coordinate of the coefficient from the vector would 
not be necessary. 

[0297] The following step 425 consists of a test aiming 
to determine whether the last coefficient of the block in 
course of being processed has been decoded. 
[0298] If so, the central processing unit goes on to 
step 435 which will determine whether the last block has 
been processed. 

[0299] If not, the central processing unit will go on to 
the following coefficient at step 430 and steps 415, 420, 
425, and 430 will be reiterated so long as the last coef- 
ficient of the block has not been processed. 
[0300] If the test at step 435 is negative, at step 440 
the central processing unit will read the following block 
in the same manner as at step 401 and reiterate the 
steps 405 to 435. 

[0301] If the last block has been processed, at step 
435 the central processing unit will go on to the calcu- 
lation of the transform that is inverse to that executed at 
the encoding, for example an inverse transform of the 
decomposition into sub-bands. 
[0302] Once this step has been accomplished, the 
method of decoding is terminated for the Image to be 
decoded. 

[0303] According to a variantform, the coordinates (p, 
q) of the first sample are read in the file. The coefficient 
of the coordinates (p, q) is then set to the value fAoptO) 
(taking into account the sign). Next the coordinates of 
the second sample are read, and at these coordinates 
the sample is set to the value fA opt (2) (taking account of 
its sign), and so on as far as the last sample encoded 
for this block. 

[0304] According to another variant form, the block 
coefficients which have not been allocated a value at 
the previous steps are allocated arbitrary values chosen 
by the decoder. In the preferred embodiment, these val- 
ues will be set at zero. 

[0305] According to Figure 15, the method of decom- 
position into sub-bands which may be contained in cir- 
cuit 51 of Figure 1 2 comprises three successive analysis 
units for decomposing the image IM into sub-bands ac- 
cording to three levels of resolution. 
[0306] In general terms, the resolution of a signal is 
the number of samples per unit length used for repre- 
senting this signal. 

[0307] The first analysis unit receives the digital im- 
age signal and applies it to two digital filters 501 and 502 
which are respectively low-pass and high -pass which fil- 
ter the image signal in a first direction, for example hor- 
izontal in the case of an image signal. After passing 
through decimators by two 510 and 520, the resulting 
filtered signals are respectively applied to two low-pass 
filters 503 and 505, and high-pass filters 504 and 506, 



which filter them in a second direction, for example ver- 
tical in the case of an image signal. Each resulting fil- 
tered signal passes through a respective decimator by 
two 530, 540, 550 and 560. The first unit outputs four 

5 sub-bands LL,, LHj, HL 1 and HH 1 of the highest reso- 
lution RES 1 In the decomposition. 
[0308] The sub-band LL, comprises the components, 
or coefficients, of low frequency, in both directions, of 
the image signal. The sub-band LH^ comprises the com- 

10 ponents of low frequency in a first direction and of high 
frequency in a second direction, of the image signal. The 
sub-band HL 1 comprises the components of high fre- 
quency in the first direction and the components of low 
frequency in the second direction. Finally, the sub-band 

to HH 1 comprises the components of high frequency in 
both directions. 

[0309] Each sub-band is an image constructed from 
the original image, which contains information corre- 
sponding to a respectively vertical, horizontal and diag- 

20 onal orientation of the image, in a given frequency band. 
[0310] The sub-band LL^ is analyzed by an analysis 
unit similar to the previous one in order to supply four 
sub-band signals LL 2 , LH 2 , HLg and HH 2 of intermediate 
resolution level RES 2 in the decomposition. The sub- 

25 band LL 2 comprises the low frequency components in 
the two directions of analysis, and Is analyzed in its turn 
by the third analysis unit which is similar to the two pre- 
vious ones. The third analysis unit supplies the sub- 
bands LL 3 , LH 3 , Hl 3 and HH 3 , of the lowest resolution 

30 level RES 3 in the decomposition, resulting from the di- 
vision into sub-bands of the sub-band LL 2 . 
[0311] Each of the sub-bands of resolution RES 2 and 
RES 3 also corresponds to an orientation in the image. 
[0312] The decomposition carried out is such that a 

35 sub-band of given resolution is divided into four sub- 
bands of lower resolution and thus has four times more 
coefficients than each of the sub-bands of lower resolu- 
tion. 

[0313] A digital image IM is represented in a schemat- 
ic ic manner in Figure 16a, whereas Figure 16b repre- 
sents the image IMD resulting from the decomposition 
of the image IM. into ten sub-bands according to three 
resolution levels. The image IMD comprises as much 
information as the original image IM, but the information 
45 is divided in frequency according to three resolution lev- 
els. 

[0314] The lowest resolution level RE S 3 comprises 
the sub-bands LL 3 , HL 3 , LH 3 and HH 3 , that is to say the 
sub-bands of low frequency in the two directions of anal- 
so ysis. The second resolution level RES 2 comprises the 
sub-bands HL 2 , LH 2 and HH 2 and the highest resolution 
level RES 1 comprises the sub-bands HL 1t LH 1 and Hh^ 
of highest frequency. 

[0315] The sub-band LL3 of lowest frequency is a re- 
55 duction of the original image. The other sub-bands are 
detail sub-bands. 

[031 6] Naturally, the number of resolution levels, and 
consequently of sub-bands, can be chosen differently, 
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for example 13 sub-bands and four resolution levels, for 
a bi-dimensional signal such as an image. The number 
of sub-bands per resolution level can also be different. 
The analysis and synthesis circuits are adapted to the 
dimension of the signal processed. 5 
[0317] Forfullerinformation on the transform Into sub- 
bands also called transform into wavelets reference 
may be made to the publication of S. Mallat, "A theory 
for multiresolution signal decomposition: the wavelet 
representation" or else to the public encoding standard 
JPEG2000 of which the description is available at the 
Internet address http://www.jpeg.org. 
[0318] According to a variant form, it is also possible 
to use the discrete cosine transform (DCT), the Fourier 
transform, or a non -I in ear transform such as a morpho- 
logical transform (ref: H. Heijmans et AL, "Multiresolu- 
tion signal decomposition schemes: part 2, Morpholog- 
ical wavelets" IEEE Transactions on Image Processing, 
Nov. 2000 ). The invention is efficient with or without de- 
composition or transformation. 
[0319] It should be noted that according to a variant 
of the invention, the low sub-band LL 3 could, due to its 
different nature, be encoded with a linear prediction 
(DPCM). 

[0320] The Image decomposed into sub-bands is next 
partitioned into blocks as shown in Figure 17. This par- 
tition is for example carried out by the circuit 51 of Figure 
12. Each detail sub-band is divided into blocks. In a pre- 
ferred embodiment the blocks will be square and of size 
32x32 sub-band coefficients, arranged so as to create 
a partition of each sub-band. However, other divisions 
into blocks are possible: blocks of any form such as rec- 
tangles, circles, non-contiguous groups of pixels, over- 
lapping or not, and may be of different size according to 
the decomposition level. 

[0321] Naturally the partitioning will be Identical and 
known by the encoder and the decoder. 
[0322] Once the image has been decomposed, the al- 
gorithms as described with reference to Figures 13 and 
1 4 will process the blocks independently, in a sequential 
or parallel manner. 

[0323] With reference to Figure 19, an embodiment 
of an encoding device 2 according to the invention is 
adapted to encode a digital signal with the object of com- 
pressing it. The encoding device is integrated into an 
apparatus, which is for example a digital camera, a dig- 
ital camcorder, a scanner, a printer, a photocopier, a fax 
machine, a database management system, or a com- 
puter. 

[0324] An image source 1 supplies a digital image to 
the encoding device 2, of which the operation will be de- 
tailed below. 

[0325] The device according to the invention compris- 
es: 

means 20 for decomposing the image into frequen- 
cy sub-bands, 

means 21 for quantizing the coefficients of the fre- 



quency sub-bands, using a predetermined set of 
quantization symbols, 

means 22 for counting the number of coefficients 
associated with each of the quantization symbols, 
means 23 for encoding the number of coefficients 
associated with each of the quantization symbols, 
means 24 for determining the location of each of the 
coefficients associated with each of the quantiza- 
tion symbols, and 

means 25 for encoding the locations of the coeffi- 
cients. 

[0326] The encoding device supplies a file containing 
data representing the compressed image, to transmis- 
sion means and/or memory storage means 3. These 
means are conventional and are not described here. 
[0327] The means 3 are connected to a decoding de- 
vice 4 according to the invention. 
[0328] The decoding device comprises: 

means 60 for decoding the number of coefficients 
associated with each of the quantization symbols, 
means 60 for decoding the location of each of the 
coefficients associated with each of the quantiza- 
tion symbols, 

means 61 for setting each of the coefficients to the 
value of the quantization symbol which corresponds 
to it. Quantized frequency sub-bands are thus 
formed. 

[0329] The decoding device 4 further comprises 
means 62 for dequantizing the quantized sub-bands. Fi- 
nally, means 63 for synthesizing the frequency sub- 
bands supply a decoded Image to a device 5 for ex- 
ploitating the decoded image. The device 5 comprises 
for example a viewing screen for displaying the decoded 
image. 

[0330] it should be noted that the encoding device and 
decoding device may be incorporated into the same ap- 
paratus, for example the computer 100 of Figure 9. 
[0331] Figure 20 shows an embodiment of a method 
of encoding an image, according to the invention. This 
method is implemented in the encoding device and com- 
prises steps S100 to S114. 

[0332] The method is carried out in the form of an al- 
gorithm which can be stored in whole or in part in any 
means of information storage capable of cooperating 
with the microprocessor. This storage means is reada- 
ble by a computer or by a microprocessor. The storage 
means is integrated or not into the device, and may be 
removable. For example, it may comprise a magnetic 
tape, a diskette or a CD-ROM (fixed memory compact 
disk). 

[0333] Step S100 is a decomposition into frequency 
sub-bands of a digital image IM to be processed accord- 
ing to the Invention. 

[0334] The signal is decomposed into frequency sub- 
bands at several resolution levels by a discrete wavelet 



15 



20 



25 



30 



35 



40 



45 



50 



17 



33 



EP 1 263 236 A2 



34 



transformation known as DWT (Discrete Wavelet Trans- 
form). Other types of decomposition can be used, for 
example Gauss-Laplace pyramids or discrete cosine 
transformations known as DCT. Hereinafter, we will con- 
sidermore particularly the case of the transformation in- 
to discrete wavelets at M resolution levels. 
[0335] A circuit for decomposition into frequency sub- 
bands, as well as an image processed by this circuit 
have been described with reference to Figures 15, 16a 
and 16b. 

[0336] The following step S1 02 is the selection of a 
first sub-band resulting from the preceding decomposi- 
tion. The sub-bands are considered one by one. 
[0337] The following step S103 is a scalar quantiza- 
tion of the current sub-band. This is conventional and 
consists overall of transforming each real conflict value 
of the sub-band into a quantization symbol. 
[0338] In a conventional manner, a predetermined 
quantization symbol set is used. A quantization symbol 
is associated with an interval of the set of real numbers 
and any coefficient belonging to this interval is quantized 
by this symbol. 

[0339] For this, a quantization step size is chosen and 
each coefficient is divided by this quantization step size 
and the result is rounded to the nearest integer. This 
integer is the quantization symbol of the coefficient con- 
sidered. 

[0340] The result of step S103 is a quantized sub- 
band, containing a set of quantization symbols respec- 
tively associated with the coefficient of the current sub- 
band. 

[0341] The following step S104 is the selection of a 
quantization symbol of the quantized sub-band. 
[0342] The following step S1 05 is a grouping together 
of all the coefficients of the sub-band associated with 
the current quantization symbol. 
[0343] At the following step S106, the number N of 
coefficients of the current group is determined. This 
number is encoded and transmitted or stored in memo- 
ry. 

[0344] The following step S1 07 is the selection of an 
initial coefficient in the group of coefficients of the sub- 
band associated with the current quantization symbol. 
[0345] At the following step S1 08, the coordinates of 
the initial coefficient are transmitted or stored in memo- 
ry. These coordinates form part of the encoding data of 
the current sub-band. For example, if the size of the cur- 
rent sub-band is.32 x 32 coefficients, a binary code of 
five bits is used to represent the x coordinate and a bi- 
nary code of five bits is used to represent the y coordi- 
nate of the initial coefficient. 

[0346] The following step S109 is the calculation of 
the shortest path for going through all the coefficients of 
the current group. This calculation is carried out accord- 
ing to a well-known technique of solving the "traveling 
salesman problem". Such techniques are for example 
described by R. E. Burkard et Al. in " Weil -solvable spe- 
cial cases of the TSP : a survey SPEZIALFORSC- 



HUNGSBEREICH F003, Berlchtn°52, December 1995. 
[0347] An example of solving this problem is also giv- 
en at the Internet address http://itp.nat.uni-magdeburg. 
de/~mertens/TSP/node2.htmi . 
5 [0348] This method is called "method by insertion" 
and relies on the principle idea of first constructing a 
path comprising a subset of the coefficients to be gone 
through. 

[0349] For example, three coefficients are chosen at 

10 random and are connected. 

[0350] Next, for each coefficient to be inserted in the 
path, calculation is made of which successive pairs of 
coefficients of the path it can be inserted so as to as to 
increase the length of the path the least possible. The 

15 cost of the coefficient under consideration is thus deter- 
mined. 

[0351] The coefficient of minimum cost is then insert- 
ed in the path. 

[0352] This procedure is reiterated until all the coeffi- 
20 cients have been inserted. 

[0353] The path so obtained is of short length, without 
it necessarily being the shortest. 
[0354] The result of step S1 09 is a set of displacement 
vectors corresponding to the path of minimum distance 
25 of the coefficients of the current group. 

[0355] The following step S110 is an entropic encod- 
ing of the displacement vectors already determined. The 
encoding used here is for example a Huffman encoding 
or an arithmetic encoding. The encoded vectors are 
30 transmitted or stored in memory. 

[0356] Thus, as a result of steps S106 to S110, the 
encoded form of a group of a sub-band comprises the 
number of coefficients of the group, the coordinates of 
the initial coefficient and the displacement vectors. As 
35 has been seen, the number of coefficients and the dis- 
placement vectors are encoded. 
[0357] As a variant, steps S1 08 to S11 0 are modified 
in the following manner. The displacement vectors of the 
coefficients are not considered, but instead the coordi- 
40 nates of each of the coefficients. These coordinates are 
then encoded according to an entropic encoding then 
transmitted or stored in memory. 
[0358] In this variant, the encoded form of a group of 
a sub-band comprises the number of coefficients of the 
45 group and the coordinates of all the coefficients of the 
group. 

[0359] The following step S1 1 1 is a test to determine 
if ail the quantization symbols used in the current sub- 
band have been considered. 

so [0360] If the response is negative, then step S111 is 
followed by the step S112 at which a following quanti- 
zation symbol is considered. Step S112 is followed by 
the previously described step S105. 
[0361] When all the symbols of the current sub-band 

55 have been processed, the response is positive at step 
S111 and this step is followed by step S113 which is a 
test to determine whether all the sub-bands have been 
processed. 
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[0362] When the response is negative, this step is fol- 
lowed by the step S114 at which a following sub-band 
is selected. Step S114 is followed by the previously de- 
scribed step S103. 

[0363] When at step S113 the response is positive, 
this means that all the sub-bands have been processed 
and the encoding of the image Is terminated. 
[0364] Figure 21 represents an embodiment of a 
method of decoding data previously encoded according 
to the method of Figure 20. 

[0365] This method is implemented in the decoding 
device and comprises the steps S221 to S233. 
[0366] The method is carried out In the form of an al- 
gorithm which can be stored in whole or in part in any 
means of information storage capable of cooperating 
with the microprocessor. This storage means Is reada- 
ble by a computer or by a microprocessor. The storage 
means is integrated or not into the device, and may be 
removable. For example, it may comprise a magnetic 
tape, a diskette or a CD-ROM (fixed memory compact 
disk). 

[0367] The data to be decoded are the encoded data 
of groups formed in the sub-bands of the Initial image. 
[0368] Step S221 is the selection of a first sub-band. 
The sub-bands are considered one by one. 
[0369] The following step S222 is the selection of a 
quantization symbol. Ail the possible quantization sym- 
bols are considered successively. 
[0370] The following step S223 is the decoding of the 
number N of coefficients of the current sub-band which 
belong to the group of the current quantization symbol. 
It should be recalled that this number forms part of the 
encoding data of the current group. The number N here 
makes it possible to know how many coefficients there 
are to decode. 

[0371] In what follows, in order to simplify, it will be 
considered that the number N is not nil, i.e. that the pos- 
sible quantization sub-bands which are considered ac- 
tually correspond to a group of the current sub-band of 
the image. 

[0372] The following step S224 is the reading of the 
coordinates of the initial coefficient of the current group. 
[0373] The following step S225 is the entropic decod- 
ing of the displacement vectors of the current group of 
coefficients. The decoding corresponds to the encoding 
which was carried out at step S1 09. 
[0374] The following step S226 is the calculation of 
the locations of the N coefficients of the current group. 
This calculation corresponds to step S108 already de- 
scribed. Based on the initial coefficient, the locations of 
the coefficients of the current group are determined. 
[0375] As a variant, when all the locations of the co- 
efficients are represented directly by their coordinates, 
then these coordinates are decoded, which determines 
the locations. 

[0376] At the following step S227, the coefficients of 
the current group are set to the value of the quantization 
symbol corresponding to them. 



[0377] The following step S228 is a test to determine 
whether the current quantization symbol is the last to be 
processed, if the response is negative, this step is fol- 
lowed by the step S229 at which a following quantization 
5 symbol is selected, so as to consider another group of 
coefficients. Step S229 is followed by the step S223 al- 
ready described. 

[0378] When the response is positive at step S228, 
this step is followed by step S230 at which the current 

w sub-band is de-quantized. 

[0379] The following step S231 is a test to determine 
whether the current sub-band is the last to be proc- 
essed. If the response is negative, this step is followed 
by the step S232 at which a following sub-band is con- 

15 sidered. Step S232 is followed by the step S222 already 
described. 

[0380] When the response is positive at step S231 , 
this step is followed by step S233 which is a reverse 
discrete wavelet transformation of the sub-bands of the 
20 image. The result is a decoded image. 



Claims 

25 1 . Method of processing a digital signal including sam- 
ples (a 1f a 2 . a 9 ) each represented by a coordi- 
nate on each axis of a reference frame with n di- 
mensions, where n is the dimension of said digital 
signal, characterized In that it includes the follow- 
so ing steps: 

determining, amongst the coordinates of said 
samples, extreme coordinates (L, R, B, T) ,on 
the x and y axes respectively of the coordinate 
35 system. 

for all possible combinations between said 
samples taken in pairs, determining a minimum 
distance amongst all the distances calculated 
along each axis between, on the one hand, the 
40 coordinate of a first sample and, on the other 

hand, the coordinate of each of the elements of 
the group comprising a second sample and 
copies of this second sample obtained by trans- 
lation of said second sample in different direc- 
ts tions and according to a value which depends 
on the extreme coordinates previously deter- 
mined. 

2. Method accordingto Claim 1 , characterized in that 
so it Includes a step of determining a minimum path 

(a 1f 84, a 6; a 7 , a 8 , a 5 , a^ a 3 , a 9 ) between the sam- 
ples from the different minimum distances previous- 
ly determined along each axis and for all the possi- 
ble combinations of pairs of samples. 

55 

3. Method accordingto Claim 2, characterized In that 
the minimum path between the samples is deter- 
mined using one of the solutions to the traveling 
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salesman mathematical problem. 

4. Method according to one of Claims 1 to 3, charac- 
terized In that the value which depends on the ex- 
treme coordinates corresponds substantially to the 
difference (L-R, R-L, B-T, T-B) between the maxi- 
mum coordinate (L t B) and the minimum coordinate 
(R f T) of the samples on each axis. 

5. Method according to one of Claims 1 to 4, charac- 
terized In that the copies of the second sample are 
made in both possible directions of each axis and 
each diagonal. 

6. Method according to one of Claims 1 to 5, charac- 
terized In that it includes a step of encoding the 
extreme coordinates (L, R, B, T) of the samples. 

7. Method according to Claim 6, characterized In that 

it includes a step of transmitting the encoded ex- 
treme coordinates (I, R, B, T) of the samples. 

8. Method according to Claim 2, characterized In that 
it includes a step of encoding the displacement vec- 
tors (W) obtained from the consecutive samples of 
the minimum path. 

9. Method according to Claim 8, characterized In 
that, for a signal of dimension 2, the coordinates 
wx, wy of a displacement vector W are supplied by 
the following formulae: 

wx = min (lx 2 - x 1 1, bc 2 + L - R - x 1 1, lx 2 - L + R - x 1 1), 

wy=min (fy^y^, lyg + B-T-yjl, ly 2 - B + T-y^), 

where x 1( y 1 and x 2 , y 2 are respectively the 
coordinates of two successive samples of a path. 

10. Method according to Claim 8 or 9, characterized 
in that it includes a step of transmitting the encoded 
displacement vectors (W) and the order of these 
vectors in the path between the samples. 

11. Method according to one of Claims 1 to 10, char- 
acterized In that the digital signal is an image sig- 
nal. 

12. Method of processing an encoded digital signal in- 
cluding a step of receiving this encoded digital sig- 
nal, characterized In that said encoded digital sig- 
nal includes, in encoded form: 

a coordinate of a so-called initial one of the 
samples of said non-encoded signal which was 
obtained along each axis of a coordinate sys- 



tem with n dimensions, where n is the dimen- 
sion of said signal, 

the extreme coordinates (L, R, B, T) on each 
axis of the samples of said non-encoded signal, 
5 - displacement vectors (W) which have each 
been obtained before coding of said signal from 
consecutive samples of a minimum path deter- 
mined between the samples of the signal, as 
well as the order of these vectors in the path. 

10 

13. Method according to Claim 12, characterized In 
that it includes the following steps of decoding: 

- the coordinate of the initial sample, 

15 - the extreme coordinates of the samples, 

- the displacement vectors (W). 

14. Method according to Claim 13, characterized In 
that it includes a step of reconstitution of the coor- 

20 dinates of each of the samples of the path, along 
each axis, according to the displacement vectors 
(W), the extreme coordinates of the samples and 
the coordinate of the initial sample. 

25 15. MethodaccordingtoClaim13or14,characterized 
In that it includes, after decoding, a step of produc- 
ing a list of the displacement vectors (W) according 
to the order of these vectors in the path. 

30 16. Method according to Claim 14, characterized in 
that it includes a step of determining the position of 
the reconstituted coordinates along each axis of 
each of the samples of the path with respect to the 
extreme coordinates of the samples. 

35 

17. Method according to Claim 16, characterized in 
that, when the determination step reveals that the 
position of a coordinate of a sample along an axis 
is situated beyond the extreme coordinates on said 

40 axis, said method includes a step of adjusting the 
position of said coordinate. 

18. Method according to Claim 17, characterized in 
that the adjustment step consists of carrying out a 

45 translation of the position of the coordinate along 
the axis in question and according to a value which 
depends on the extreme coordinates on said axis 
and the position of the coordinate with respect to 
said extreme coordinates. 

50 

19. Method according to Claim 18, characterized in 
that the value which depends on the extreme coor- 
dinates on the axis corresponds substantially to the 
difference (L-R, R-L, B-T, T-B) between the maxi- 

55 mum coordinate (L, B) and the minimum coordinate 
(R, T) of samples on the axis in question. 

20. Method according to one of Claims 17 to 19, char- 
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acterlzed In that, for a signal of dimension 2. the 
step of adjusting the position of the coordinates x 2 
and y 2 obtained from the coordinates Xj and y 1 of 
the previous sample on the path is given by the fol- 
lowing formulae: 

If x 1 + wx < L, x 2 » x 1 + wx + R - L 



If x t + wx > R, x 2 = x 1 + wx - R + L, 



lfy 1 +wy>T, y 2 =:y 1 +wy + B-T, 



If y 1 + wy < B, y 2 = y 1 + wy - B + T, 

where wx and wy are the coordinates of the dis- 
placement vector W in question. 

21 . Device for processing a digital signal including sam- 
ples (a 1: a 2 , Qq) each represented by a coordi- 
nate on each axis of a coordinate system with n di- 
mensions, where n is the dimension of said digital 
signal, characterized in that it includes: 

means of determining, amongst the coordi- 
nates of said samples, extreme coordinates (L, 
R, B, T) on the x and y axes respectively of the 
coordinate system. 

for all possible combinations between said 
samples taken in pairs, means of determining 
a minimum distance amongst ail the distances 
calculated along each axis between, on the one 
hand, the coordinate of a first sample and, on 
the other hand, the coordinate of each of the 
elements of the group comprising a second 
sample and copies of this second sample ob- 
tained by translation of said second sample in 
different directions and according to a value 
which depends on the extreme coordinates 
previously determined. 

22. Device according to Claim 21, characterized in 
that it has means of determining a minimum path 
(a 1( a 4 , ae, a 7 , a 8> a 5 , a 2 , a 3 , ag) between the sam- 
ples from the different minimum distances previous- 
ly determined along each axis and for all possible 
combinations of pairs of samples. 

23. Device according to Claim 22, characterized in 
that the means of determining a minimum path be- 
tween the samples use one of the solutions to the 
traveling salesman mathematical problem. 

24. Device according to one of Claims 21 to 23, char- 
acterized In that the value which depends on the 



25. 



extreme coordinates corresponds substantially to 
the difference (L-R, R-L, B-T, T-B) between the 
maximum coordinate (L, B) and the minimum coor- 
dinate (R, T) of the samples on each axis. 

Device according to one of Claims 21 to 24, char- 
acterized In that it includes means of encoding the 
extreme coordinates (L, R), (B, T) of the samples. 



10 26. Device according to Claim 25, characterized In 
that It includes means of transmitting the extreme 
coordinates (L, R), (B, T) of the encoded samples. 

27. Device according to Claim 22, characterized in 
15 that it has means of encoding the displacement 

vectors (W) obtained from the consecutive samples 
of the minimum path. 

28. Device according to Claim 27, characterized in 
20 that it has means of transmitting the encoded dis- 
placement vectors (W) and the order of these vec- 
tors In the path between the samples. 

29. Device according to one of Claims 21 to 28, char- 
ts acterlzed in that the digital signal is an image sig- 
nal. 

30. Device for processing an encoded digital signal 
having means of receiving this coded digital signal, 

30 characterized in that the encoded digital signal in- 
cludes, in encoded form: 



a coordinate of a so-called initial one of the 
samples of said non-encoded signal which was 
obtained along each axis of a coordinates sys- 
tem with n dimensions, where n is the dimen- 
sion of said signal, 

the extreme coordinates (L, R, B, T) on each 
axis of the samples of said non-encoded signal, 
displacement vectors (W) which have each 
been obtained before encoding of said signal 
from the consecutive samples of a minimum 
path determined between the samples of the 
signal, as well as the order of these vectors in 
the path. 
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31. Device according to Claim 30, characterized In 
that it includes means of decoding: 

- the coordinate of the initial sample, 
the extreme coordinates of the samples, 
the displacement vectors (W). 

32. Device according to Claim 31, characterized in 
that it includes a means of reconstitution of the co- 
ordinates of each of the samples of the path, along 
each axis, according to the displacement vectors 
(W), the extreme coordinates of the samples and 
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puter program for executing the steps of the method 
of processing a digital signal according to one of 
Claims 1 to 11. 

5 43. Partially or totally removable information storage 
means which can be read by a computer or a mi- 
croprocessor containing code instructions of a com- 
puter program for executing the steps of the method 
of processing an encoded digital signal according 
10 to one of Claims 1 2 to 20. 



41 

the coordinate of the initial sample. 

33. Device according to Claim 31 or 32, 
In that it has means of producing a list of the dis- 
placement vectors (W) according to the order of 
these vectors in the path. 

34. Device according to Claim 32, characterized In 
that it has means of determining the position of the 
reconstituted coordinates along each axis of each 
of the samples of the path with respect to the ex- 
treme coordinates of the samples. 

35. Device according to Claim 34, characterized in 
that it has means of adjusting the position of a co- 
ordinate of a sample along an axis which are used 
when the determination means reveal that the po- 
sition of the said coordinate is situated beyond the 
extreme coordinates on said axis. 

36. Device according to Claim 35, characterized In 
that the adjustment means perform a translation of 
the position of the coordinate along the axis in ques- 
tion and according to a value which depends on the 
extreme coordinates on the said axis and the posi- 
tion of the coordinate with respect to the said ex- 
treme coordinates. 

37. Device according to Claim 36, characterized in 
that the value which depends on the extreme coor- 
dinates on the axis corresponds substantially to the 
difference (L-R, R-L, B-T, T-B) between the maxi- 
mum coordinate (L, B) and the minimum coordinate 
(R, T) of the samples on the axis in question. 

38. Data processing apparatus, characterized In that 

it includes a digital signal processing device accord- 
ing to one of Claims 21 to 29. 

39. Data processing apparatus, characterized In that 
it has a device for processing an encoded digital sig- 
nal according to one of Claims 30 to 37. 

40. information storage means which can be read by a 
computer or a microprocessor containing code in- 
structions of a computer program for executing the 
steps of the digital signal processing method ac- 
cording to one of Claims 1 to 11. . _ . . 

41. Information storage means which can be read by a 
computer or a microprocessor containing code in- 
structions of a computer program for executing the 
steps of the method of processing an encoded dig- 
ital signal according to one of Claims 12 to 20. 

42. Partially or totally removable information storage 
means which can be read by a computer or a mi- 
croprocessor containing code instructions of a com- 



44. Computer program which can be loaded into a pro- 
grammable apparatus, characterized In that it 
contains sequences of Instructions or portions of 

15 software code for implementing the steps of the 
method of processing a digital signal according to 
one of Claims 1 to 11 , when this computer program 
is loaded and executed by the programmable ap- 
paratus. 

20 

45. Computer program which can be loaded into a pro- 
grammable apparatus, characterized In that it 
contains sequences of instructions or portions of 
software code for implementing the steps of the 

25 method of processing an encoded digital signal ac- 
cording to one of Claims 12 to 20, when this com- 
puter program is loaded and executed by the pro- 
grammable apparatus. 

30 46. Method of encoding a digital signal composed of a 
set of samples, with each sample there being asso- 
ciated an amplitude and at least one coordinate in 
a coordinate system, the method being character- 
ized in that it comprises the steps of: 

35 

classifying at least some of the samples ac- 
cording to a predetermined criterion, 
obtaining a curve representing the variation in 
the samples classified according to the prede- 
40 termined criterion, 

encoding the data representing the curve ob- 
tained, 

encoding information representing at least one 
coordinate of at least one coefficient in an order 
45 obtained from the classification. 

47. Method according to Claim 46, characterized in 

that the digital signal consists of at least one digital 
image. 

50 

48. Method according to Claim 47, characterized in 

that the digital image is first divided into sub-parts. 

49. Method according to Claim 48, characterized in 
55 that the division into sub-parts is carried out accord- 
ing to a decomposition into sub-bands. 

50. Method according to Claim 49, characterized In 
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that the sub-parts are the detail sub-bands, the sub- 
bands being divided into blocks of predetermined 
size and the classification of at least some of the 
samples being carried out from among the samples 
of each block. 

51. Method according to one of Claims 46 to 50, char- 
acterized in that the predetermined criterion is the 
amplitude of the samples. 

52. Method according to Claim 51, characterized in 
that the method further comprises an encoding of 
the sign of the amplitude of the coefficients, and the 
predetermined criterion is the absolute value of the 
amplitude of the samples. 

53. Method according to Claim 52, characterized In 
that only the samples whose amplitude is greater 
than a predetermined threshold will be classified. 

54. Method according to any one of Claims 46 to 53, 
characterized in that the curve representing the 
variation according to the predetermined criterion 
of the samples classified is chosen from among a 
predetermined family. 

55. Method according to Claim 54, characterized in 
that the predetermined family is the family of de- 
creasing exponentials, the choice of the curve 
among the family of exponentials being made ac- 
cording to the minimization of the root mean square 
error. 

56. Method according to any one of Claims 46 to 55, 
characterized in that the information representing 
at least one coordinate of at least one coefficient 
are the coordinates of the coefficients in the image. 

57. Method according to one of Claims 46 to 56, char- 
acterized in that the information representing at 
least one coordinate of at least one coefficient are 
the vectors of displacement between two coeffi- 
cients which are close according to the classifica- 
tion. 

58. Method according to Claim 57, characterized in 
that the displacement vectors are encoded accord- 
ing to a Huffman code, constraints of rate of the en- 
coded signal and/or of distortion of the decoded dig- 
ital signal being taken into account in determining 
a travel path of the samples. 

59. Method of transmitting a digital signal which has 
been encoded according to an encoding method ac- 
cording to Claims 46 to 58, characterized In that 
the encoded data representing the curve obtained 
and the encoded information representing at least 
one coordinate are transmitted via a network. 



44 

60. Method of decoding a digital signal comprising co- 
efficients, the method being characterized in that 
ft comprises the steps of: 

5 - decoding information representing a curve of 
variation of the amplitude of the classified co- 
efficients, 

decoding information representing at least one 
coordinate of the coefficients in a predeter- 
10 mined order and associating with the coeffi- 

cients an amplitude obtained from the informa- 
tion representing the decoded curve and the or- 
der of the coefficients. 

15 61. Method according to Claim 60, characterized in 
that the digital signal is constituted by at least one 
digital image. 

62. Method according to Claim 61, characterized in 
20 that it further comprises a step of transformation 

that is Inverse to decomposition into sub-bands of 
the image. 

63. Method according to Claim 62, characterized in 
25 that the sub-parts are the detail sub-bands, the de- 
tail sub-bands being divided into blocks of prede- 
termined size and the decoding taking place block 
by block. 

30 64. Method according to any one of Claims 60 to 63, 
characterized in that the information representing 
at least one coordinate of at least one coefficient 
are the spatial coordinates of coefficients in the im- 
age. 

35 

65. Method according to one of Claims 60 to 63, char- 
acterized In that the information representing at 
least one coordinate of at least one coefficient are 
the vectors of displacement between two coeffi- 

40 cients which are close according to the classifica- 
tion. 

66. Method according to Claim 60, characterized in 
that it further comprises a step of decoding of the 

45 sign of the amplitude of the coefficients and In that 
the curve of variation is a curve of variation of the 
absolute value of the amplitude of the coefficients 
classified. 

so 67. Device for encoding a digital signal composed of a 
set of samples, with each sample there being asso- 
ciated an amplitude and at least one coordinate in 
a coordinate system the device being character- 
ized in that it comprises: 

55 

means of classifying at least some of the sam- 
ples according to a predetermined criterion, 
- means of obtaining a curve representing the 
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variation according to the predetermined crite- 
rion in the samples classified, 
means of encoding data representing the curve 
obtained, 

means of encoding information representing at 
least one coordinate of at least one coefficient 
in an order obtained from the classification. 

68. Device according to Claim 67, characterized in 
that the digital signal is constituted by at least one 
digital image. 

69. Device according to Claim 68, characterized in 
that it further comprises means of dividing the im- 
age into sub-parts. 

70. Device according to Claim 69, characterized in 
that the means of dividing the image are means of 
decomposition into sub-bands. 

71. Device according to Claim 70, characterized in 
that the sub-parts are the detail sub-bands, the sub- 
bands being divided into blocks of predetermined 
size and the classification of at least some of the 
samples being carried out among the samples of 
each block. 

72. Device according to one of Claims 67 to 71, char- 
acterized in that the predetermined criterion is the 
amplitude of the samples, 

73. Device according to Claim 72, characterized in 
that the method further comprises an encoding of 
the sign of the amplitude of the coefficients, and the 
predetermined criterion is the absolute value of the 
amplitude of the samples. 

74. Device according to Claim 72, characterized in 
that only the samples whose amplitude is greater 
than a predetermined threshold will be classified. 

75. Device according to any one of Claims 72 to 74, 
characterized in that the curve representing the 
variation according to the predetermined criterion 
of the samples classified is chosen from among a 
predetermined family. 

76. Device according to Claim 67, characterized in 
that the predetermined family is the family of de- 
creasing exponentials, the choice of the curve 
among the family of exponentials being carried out 
by means of selections as a function of criteria of 
minimization of the root mean square error. 

77. Device according to any one of Claims 67 to 76, 
characterized in that the information representing 
at least one coordinate of at least one coefficient 
are the coordinates of the coefficients in the image. 



78. Device according to one of Claims 67 to 77, char- 
acterized in that the information representing at 
least one coordinate of at least one coefficient are 
the vectors of displacement between two coeffi- 

5 cients which are close according to the classifica- 
tion. 

79. Device according to Claim 76, characterized in 
that the displacement vectors are encoded accord- 

10 ing to a Huffman code, a travel path of the samples 
being determined taking into account constraints of 
rate of the encoded signal and/or of distortion of the 
decoded digital signal. 

*5 80. Device for transmitting a digital signal which has 
been encoded according to an encoding method ac- 
cording to Claims 67 to 79, 
characterized in that the encoded data represent- 
ing the curve obtained and the encoded information 
20 representing at least one coordinate are transmitted 
via a network. 

81 . Device for decoding a digital signal comprising co- 
efficients, the device being characterized in that it 

25 comprises: 

means of decoding information representing a 
curve of variation of the amplitude of the coef- 
ficients classified, 
30 - means of decoding information representing at 
least one coordinate of the coefficients in a pre- 
determined order and of associating with the 
coefficients an amplitude obtained from infor- 
mation representing the decoded curve and the 
35 order of the coefficients. 

82. Device according to Claim 81, characterized in 
that the digital signal is constituted by at least one 
digital image. 

40 

83. Device according to Claim 82, characterized in 
that It further comprises means of transformation 
that is inverse to decomposition into sub-bands of 
the image. 

45 

84. Device according to Claim 83, characterized in 
that the sub-parts are the detail sub-bands, the de- 
tail sub-bands being divided into blocks of prede- 
termined size and the decoding takes place block 

so by block. 

85. Device according to any one of Claims 82 to 84, 
characterized In that the information representing 
at least one coordinate of at least one coefficient 

55 are the spatial coordinates of coefficients in the im- 
age. 

86. Device according to one of Claims 81 to 84, char- 
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35 



40 



45 



24 
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acterlzed In that the information representing at 
least one coordinate of at least one coefficient are 
the vectors of displacement between two coeffi- 
cients which are close according to the classifica- 
tion. 

87. Device according to Claim 81 , characterized In 
that it further comprises means for decoding the 
sign of the amplitude of the coefficients and In that 
it is adapted to consider a curve of variation which 
is a curve of variation of the absolute value of the 
amplitude of the coefficients classified. 

88. Storage medium storing a program for implement- 
ing the method according to any one of Claims 46 
to 66. 

89. Storage medium according to claim 88, character- 
ised in that said storage medium is detachably 
mountable on a device according to any one of 
Claims 67 to 87. 

90. Storage medium according to claim 88 or 89, char- 
acterised in that said storage medium is a floppy 
disk or a CD-ROM. 

91 . Computer program on a storage medium and com- 
prising computer executable instructions for caus- 
ing a computer to encode a digital signal according 
to any one of claims 46 to 58. 

92. Computer program on a storage medium and com- 
prising computer executable instructions for caus- 
ing a computer to transmit a digital signal according 
to claim 59. 

93. Computer program on a storage medium and com- 
prising computer executable instructions for caus- 
ing a computer to decode a digital signal according 
to any one of claims 60 to 66. 

94. A method of encoding a set of data representing 
physical quantities, the set of data comprising coef- 
ficients, 

characterized in that it comprises the steps 

of: 

quantizing (S103) the coefficients using a pre- 
determined set of quantization symbols, 
counting (S106) the number of coefficients as- 
sociated with each of the quantization symbols, 
encoding (S1 06) the number of coefficients as- 
sociated with each of the quantization symbols, 
determining (S108, S109) the location of each 
of the coefficients associated with each of the 
quantization symbols, and 
encoding (S110) the locations of the coeffi- 
cients. 



48 

95. A method of encoding according to Claim 94, char- 
acterized In that the locations of the coefficients 
are determined by: 

5 - the coordinates of an initial coefficient (S1 08) 
in the set of data, 

displacement vectors (S109) for the other co- 
efficients, the displacement vector of a coeffi- 
cient being calculated from another coefficient 
10 associated with the same symbol. 

96. A method of encoding according to Claim 95, char- 
acterized In that the encoding of the locations of 
the coefficients comprises (S110) an entropic en- 

15 coding of the displacement vectors, 

97. A method of encoding according to Claim 94, char- 
acterized In that the locations of the coefficients 
are represented by their coordinates in the set of 

20 data. 

98. A method of encoding according to Claim 97, char- 
acterized In that the encoding of the locations of 
the coefficients comprises an entropic encoding of 

25 the coordinates of the coefficients. 

99. A method of encoding according to any one of 
Claims 94 to 98, characterized in that the set of 
data is a frequency sub-band signal resulting from 

30 a decomposition (S1 01 ) into sub-bands of an initial 
signal. 

100. A method of decoding a set of data representing 
physical quantities encoded by the method accord- 

35 ing to any one of Claims 94 to 99, 

characterized in that it comprises the steps 

of: 

decoding (S223) the number of coefficients as- 
40 sociated with each of the quantization symbols, 

- decoding (S224, S225, S226) the location of 
each of the coefficients associated with each of 
the quantization symbols, 
setting (S227) each of the coefficients to the 
45 value of the quantization symbol which corre- 

sponds to it. 

101 .A method of decoding according to Claim 100, 
characterized in that the decoding of the location 
so of a coefficient comprises (S225) the decoding of 
the displacement vector corresponding to it. 

102. A method of decoding according to Claim 100, 
characterized in that the decoding of the location 

55 of a coefficient comprises the decoding of its coor- 
dinates. 

103. A device for encoding a set of data representing 
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physical quantities, the set of data comprising coef- 
ficients, 

characterized In that it comprises: 

means (21) for quantizing the coefficients using 5 
a predetermined set of quantization symbols, 
means (22) for counting the number of coeffi- 



cients associated with each of the quantization 
symbols, 

- means (23) for encoding the number of coeffi- 
cients associated with each of the quantization 
symbols, 

means (24) for determining the location of each 
of the coefficients associated with each of the 
quantization symbols, and 
means (25) for encoding the locations of the co- 
efficients. 

104. An encoding device according to Claim 103, char- 
acterized In that the means (24) for determining 
the locations of the coefficients are adapted to de- 
termine them by: 

the coordinates of an initial coefficient in the set 
of data, 

displacement vectors for the other coefficients, 
the displacement vector of a coefficient being 
calculated from another coefficient associated 
with the same symbol. 

105. An encoding device according to Claim 104, char- 
acterized In that the means (25) for encoding the 
locations of the coefficients are adapted to perform 
an entropic encoding of the displacement vectors. 

106. An encoding device according to Claim 103, char- 
acterized In that the means (24) for determining 
the locations of the coefficients are adapted to de- 
termine them by their coordinates in the set of data. 

107. An encoding device according to Claim 106, char- 
acterized in that the means (25) for encoding the 
locations of the coefficients are adapted to perform 
an entropic encoding of the coordinates of the co- 
efficients. 

108. An encoding device according to any one of Claims 
103 to 107, characterized in that it is adapted to 
process a set of data which is a frequency sub-band 
signal resulting from a decomposition (20) into sub- 
bands of an initial signal. 

109. A device for decoding a set of data representing 
physical quantities encoded by the device accord- 
ing to any one of Claims 1 03 to 1 08, 

characterized in that it comprises: 

means (60) for decoding the number of coeffi- 
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cients associated with each of the quantization 
symbols, 

means (60) for decoding the location of each of 
the coefficients associated with each of the 
quantization symbols, 

means (61 ) for setting each of the coefficients 
to the value of the quantization symbol which 
corresponds to it. 

10 110.A decoding device according to Claim 109, char- 
acterized In that the means (60) for decoding the 
location of a coefficient are adapted to perform the 
decoding of the displacement vector corresponding 
to it. 

15 

111. A decoding device according to Claim 109, char- 
acterized in that the means (60) for decoding the 
location of a coefficient are adapted to perform the 
decoding of its coordinates. 

20 

112. An encoding device according to any one of claims 
103 to 108, 

characterized in that the means for quantization, 
counting, determination and encoding are incorpo- 
25 rated in: 

a microprocessor (1 03), 
a read only memory (104), comprising a pro- 
gram for processing the data, and 
30 - a random access memory (106) comprising 
registers adapted to record variables modified 
during the execution of said program. 

113. A decoding device according to any one of claims 
35 109 to 111, 

characterized in that the means for decoding and 
setting to the value are incorporated in: 

a microprocessor (1 03) , 
40 - a read only memory (104), comprising a pro- 
gram for processing the data, and 
a random access memory (106) comprising 
registers adapted to record variables modified 
during the execution of said program. 

45 

1 1 4. A device for processing (1 00) a digital image, char- 
acterized In that it comprises means adapted to 
implement the method according to any one of 
Claims 94 to 102. 

50 

1 15. A device for processing (1 00) a digital image, char- 
acterized In that it comprises the device according 
to any one of Claims 103 to 113. 

55 1 16. Storage medium storing a program for implement- 
ing the method according to any one of Claims 94 
to 102. 
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117.Storage medium according to claim 116, charac- 
terised in that said storage medium is detachably 
mountable on a device according to any one of 
Claims 67 to 87. 

5 

118.Storage medium according to claim 116 or 117, 
characterised in that said storage medium is a 
floppy disk or a CD-ROM. 

1 1 9. Computer program on a storage medium and com- 10 
prising computer executable instructions for caus- 
ing a computer to encode a set of data according to 
any one of claims 94 to 99. 

1 20. Computer program on a storage medium and com- 15 
prising computer executable instructions for caus- 
ing a computer to decode a set of data according to 
any one of claims 1 00 to 1 02. 

20 
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